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1.0 Introduction and Summary 

Reliable evaluation of the dynamic characteristics of large structural 
systems such as the Space Shuttle vehicle will require combined analytical and 
experimental efforts. Because of its size and physical configuration, full-scale 
vibration testing of the entire vehicle may be impractical. Testing 
would then be limited to isolated substructures of the system such as the 
Orbiter, Tank and Booster. Analytical means are being sought to utilize data 
from sub-structure tests to predict characteristics of the connected assembly. 
Modal coupling or synthesis has been used in the past to predict undamped dynamic 
characteristics including vibration modes and frequencies. While the Shuttle 
will require special consideration in the application of modal synthesis because 
of its unusual configuration, existing methods have been developed to a point 
where the capability is within the state of the art. 

This is not true in the case of damping, however. A generally applicable 
method for predicting the damping in a structural system on the basis of substruc- 
ture damping has not yet been demonstrated. Yet, a reasonably accurate assessment 
of structural damping will be necessary in order to properly evaluate conditions 
of resonant response to either sinusoidal or random excitation. If reliable 
predictions are to be made without recourse to testing the entire vehicle, this 
capability will be essential. 

Reported herein is the development of a method for damping synthesis. 

The method has been applied successfully to ..eal structures as well as analytical 
models. It depends on the ability to determine an appropriate modal damping 
matrix for each substructure. Previous attempts at using this approach have 
failed because of invalid restrictions placed on the nature of substructure damp- 
ing. In the past, modal damping matrices have been assumed diagonal for lack of 
being able to determine the coupling terms which are significant in the general 
case of nonproportional damping. This problem has been overcome by formulating 
the damped equations of motion e? a linear perturbation of the undamped equations 
for lighc structural damping. Damped modes are defined as complex vectors 
derived from the complex frequency response vectors of each substructure and are 
obtained directly from sinusoidal vibration tests. The damped modes are used 
to compute first order approximations to the modal damping matrices. The pertur- 
bation approach avoids ever having to solve a complex eigenvalue problem. This 
is a useful fact which a number of investigators have failed to recognize. 
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2. Background 

A variety of dissipative mechanisms contribute to the overall damp 4 -' - 
properties of a complex structure. They may Include material damping cau- by 
microscopic slip between particles within the material, thermal loss due to molec- 
ular abrasion and, in most cases, structural hysteresis. At the macroscopic level, 
losses may arise from relative slip in mechanical joints, from acoustic radiation, 
and from viscous losses due to liquid sloshing or vibration of viscoelastic 

materials. Detailed investigations of damping in simple structures have revealed 

r 

various nonlinear damping laws [1,2] . Some lead to linear equations of motion, 
while others do not. However, specific nonlinear phenomena are likely to be 
obscured when appearing together in a complex structure. One may find it difficult 
to define a particular nonlinear damping law for general application. Although 
Chang [3] has concluded that a simple nonlinear damping lav may be appropriate 
for some complex structures, based on his interpretation of test data from Saturn- 
type vehicles and their scale models, other interpretations of the same data 
indicate linear dmaping with respect to amplitude for many of the modes considered 
individually. Chang's interpretation ignores differences in damping among the 
various modes, with no clear physical Justification. 

Since modal synthesis presupposes linearity, it is consistent to treat 
damping in the same manner. In view of the computational advantages offered by 
linearization and the lack of much strong physical evidence supporting the 
selection of a particular nonlinear model, the tentative assumption of equivalent 
linear damping appears to be justified. 

Recent attempts have been made to predict the modal damping properties 
of structural systems based on the modal damping of their respective substructures 
when damping is known to be linear. Collins, Hart, Hurty, and Kennedy [4] showed 
that the uncoupled modal damping properties of substructures do not constitute 
sufficient information to enable one to predict damping in the system modes using 
the standata matrix transformation approach. Kana ar»J Huzar [5] confirmed this 
finding and went on to develop a completely different approach based on total 
energy dissipation rather than energy dissipation associated with any particular 
mode. They adopted a representation of structural damping similar to Chang's, 


Numbers in square brackets d -ignate references ac end of text. 
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plotting dissipative energy per cycle versus kinetic energy. These quantities, 
obtained experimentally at each of the structure's resonant frequencies within a 
certain range, define points through which smooth curves are drawn. The curves 
are used to define the overall dissipative properties of the structure through a 
continuous frequency spectrum. A^ter determining the coupred system modes, the 
kinetic energy of each component is determined. Going back to the energy plots, 
one may pick off appropriate levels of dissipative energy for each component and 
sum them to get the total dissipative energy of the system. From this, damping 
is found for each of the system modes. While the authors demonstrate successful 
application of the method for the structural models employed, certain pitfalls are 
noted. A significant one has to do with how the points fall on the energy plots. 
If they do not align themselves reasonably well along some smooth curve, obvious 
difficulties arise. This could very well be the case when the amount of damping 
varies greatly from mode to mode, particularly if the variability is irregular. 

In this case, the authors suggest the possible grouping of points into different 
classes of modes in the hope that smooth curves may still be drawn. The chance 
for ambiguity here is not difficult to imagine, and further study will be required 
before the method is acceptable for general use. 

On the other hand, the matrix approach is not affected by nonuniform 
in the distribution of damping. In addition, it lends itself to a form common 
used by analysts. It is clear (at least in the case of linear viscous damping) 
that failure to predict damping in the system modes on the basis of component 
modal damping is attributable to the fact that the off-diagonal elements of 
the modal damping matrix are neglected. Hasselman [6] showed that the diagonal 
elements of the modal damping matrix correspond to the modal damping constants 
for each mode. Furthermore, it is shown that, in general, the off-diago..al 
elements will be of the same order of magnitude as the diagonal elements. They 
may therefore not be neglected in transformations involving the component modal 
damping matrices. The problem is to somehow determine these elements. 
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3. Objectives and Scope 

This investigation is part of a broader study directed toward the devel- 
opment of modal coupling procedures for the Space Shuttle vehicle. The principal 
objectives have been to formulate analytical procedures and then verify them using 
experimental data. TRW and the Grumman Aerospace Corporation (GAC) have been 
jointly responsible for the fomulation of procedures. Experimental verification 
has been made by TRW and GAC on the baslr. sf test data furnished by the NASA 
Langley Research Center (LaRC). GAC has performed that part of the study which 
involves synthesis of system modes and frequencies, while TRW has developed and 
Implemented a damping synthesis technique, *ihich is discussed in this report. 

The scope of the present work fits within the framework of conventional 
modal synthesis methods. Modal damping matrices derived at the substructure level 
are operated on by the seme transformations used to couple and diagonalize the 
undamped equations of motion at the system level. The main focus is placed on 
obtaining the appropriate modal damping matrices. 

Some basic assumptions are made: that the damping forces in the struc- 

ture are small compared to either the elastic forces or inertial forces, and that 
they vary linearily with respect to the amplitude of the motion. It has further- 
more been assumed in the present development that the damping in the structure may 
be represented by constant viscous damping. The Voigt model is used (as opposed 
to the Maxwell model, for example). Thi* is clearly an appropriate first step to 
take in the investigation of a linearized matrix approach. 

To relate this study to other recent work done in the same area, it is 
recalled that both [4] and [5] reported attempts to use the standard matrix coup- 
ling transformations as a means of predicting modal damping at the system level 
based on measurements of uncoupled modal damping at the substructure level. That 
is, diagonal modal damping matrices were assumed. Reference [4] concluded that, 
in general, more information is required to enable the prediction of damping in 
this way. Reference {5} concluded somewhat to the contrary that » different 
method should be sought, and proceeded to develop an energy approach instead. 

The present work begins with the conclusions of Reference (4] and movec in the 

direction of Identifying and obtaining the additional information required by 

the matrix method. The basic problem can be formulated so that within the scope 

of the linearized approach, information is provided by the off-diagonal terms of 

the substructure modal damping matrices and a description of any damping which 

may take place between substructural interfaces. ' 
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Corresponding to the major objectives of the effort, this report is 
organized into three parts: A, B, and C. Part A includes the formulation of 

methods. Part B includes the presentation and discussion of results obtained by 
applying the theory to both ideal and practical problems. Application of the 
theory to ideal problems which satisfy the inherent assumptions serves to verify 
the theoretical development and provide additional insight into practical appli- 
cation. As such, it is an essential link between the theoretical work which tends 
to become involved, and the results of practical application which embody a number 
of uncertainties associated with data acquisition and reduction. 

While every attempt has been made to minimize the Influence of these 
uncertainties, a systematic treatment of their influence on the fi'ial results was 
not within the scope of this study. Vibration tests were conducted with standard 
equipment, and data were recorded on analog tape. Analog data reduction proce- 
dures were also used. Manual transfer of the analog data to digital form on 
punched cards was employed- Although the analog data were visually displayed in 
digital form, some real-time variation did occur so that the person recording the 
data was required, in some cases, to mentally average the results before writing 
them down. 

The need for a final evaluation of the study , including a discussion 
of the current applicability of methods, suggested the addition of Part C. Some 
general comments pertaining to the matrix method and a comparison with the energy 
method of Kana and Huzar are ’made. Conclusions and recommendations are contained 
therein. 
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4. Basic Equations for Modal Synthesis 

There are various coupling procedures available for modal synthesis 
depending on the type of substructure modes used. From the standpoint of damping 
synthesis, the basic equations are the same. Without going into detail, the gen- 
eral procedure may be outlined rather simply. 

It is convenient to begin writing equations for each substructure in terms 
of a discrete coordinate system in which the displacement vector will be denoted 
as x*. These coordinates may be thought of as resulting from some finite element 
discretization of the substructure, although it is recognized that the equations 
in this form need not be defined when basing an analysis on vibration test data. 
However, a mass matrix corresponding to such a coordinate system will be required 
so that conceptually this step is necessary. The equations may be written 


p x + p 


IS. A 


( 1 ) 


The force vector f* corresponds to displacements x*, while p*, p*, and k* are 
square matrices containing, respectively, the mass, damping, and stiffness coef- 
ficients of the equations. In the case of free vibration , f* will represent 
only the interaction forces between adjacent components. Similar sets of equa- 
tions may be written for each of N components and arrayed in a diagonal pattern 
such that a complete set of equations representing the complete structure becomes 


where 


x 



pX + P X + KX * f 


and p * 



( 2 ) 


The matrices p and k are of the same form as p. 

There are three transformations involved in going from (2) to a reduced 
set of equations in system modal coordinates. The first transformation, denoted 
by ij>^, involves the hypothetical undamped substructure iodes. The subscript R is 
used to signify that these modes are real. The transformation may be written 


■■'if 


(3) 



T 


where 
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and <J> is a rectangular matrix whose columns define the characteristic deformation 

shapes of the substructure. They may include static as well as dynamic shapes. 

The second transformation is a rectangular matrix denoted by B, and is 

used to introduce compatibility constraints at substructural interfaces. The 

final transformation. 4> , involves the hypothetical undamped modes of the system 

where each column of is a system eigenvector in the coupled component mode 

R 

coordinates. These two transformations relate the p coordinate vector to a new 
coordinate vector ri, whose elements correspond to tne system eigenvectors ♦ • 


P - B$ R n * Tn. 


It is convenient both conceptionally and notationally to combine B and $ R into 
the single transformation matrix T. Transformation of (2) to the n coordinate 
system gives equations of the form 


In + Cn + An * 0 


when $ is normalised to give unit modal mass. In other words, 
K 

T T 

t * R u * r t - i 

T T 

T \ P *R T * C 


T T 

T * R k* r T - A 


o 


where I is an identity matrix, A is a diagonal matrix of system eigenvalues 
(frequencies squared), and C is, in general, a fully populated modal damping 
matrix at the system level. 
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The nonhomogeneous form of Equation (5) is used to solve the forced 
o response problem. The off-diagonal elements of C are usually neglected while the 

diagonal elements, according to [6] are given by 

1/2 


CjJ - 2 Cj Aj 


( 7 ) 


; being the critical damping ratio for the jth mode. Without full-scale testing 
of the composite structure, cannot be obtained directly. The problem then is 
to find some other means of determining C^j. 

T 

In Equation (6b) , it is convenient to let = c where 

1 


’ N 

L c 

i i^ i i 

and c = K P 4> d * The scalar C.. may then be expressed in the form 
R R j j 


'jj 


E t j 


T 

i i~i 
c T. 


( 8 ) 


where T^ denotes the jth column of T and 


In , 

T J 

corresponding to the partitioning of c. The matrices T* are given by synthesis 
of the undamped equations. It remains then to #valuate c^. 


i 
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5. Forms of the Substructure Modal Damping Matrix 

Identification of substructure damping properties will depend on the 
type of substructure modes used in the synthesis. The mode types are governed 
by the boundary conditions imposed. Interface boundaries may be free, fixed, 
or reflect some intermediate degree of fixity by (by application of mass loading, 
for example). The case involving free-interface boundary conditions is the sim- 
plest , both conceptually and from an experimental viewpoint. In this case, the 

modal matrix will contain rigid body and elastic modes. For free vibration 
i R 

modes, $ R may be partitioned to distinguish between the rigid body and elastic 
(normal) modes. 




i i 

Relative to the modal coordinates p , c may be partitioned as 


■‘i [c : c] 


iRR , iRN 
c 1 c 

' Inr "« "iNN' 


Iu the absence of external damping. 


The submatrix c may be determined by the method discussed In [6]. After 
obtaining til the c* matrices, i ■ 1 to N, modal damping for the coupled system 
can be determined by Equation (8). For mass-loaded interface modes, the form of 

± 

c is . 1st given by (10) . 

When fixed- interface substructure modes are used, the determination of 
c~ will be more difficult. In this case, $ may include three different types 
of modes — rigid body, constraint, and normal modes. The constraint modes appear 
oniy for redundantly interconnected substructures and are introduced so that the 
motion ot interface boundary points is completely defined. These are usually 
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chosen to be static deformation shapes which are obtained by displacing each 
boundary coordinate sequentially while holding all others fixed and applying 
no loading whatsoever at points other than attachment points. Under these condi- 
tions, <J>* may be partitioned in the manner 


where <J> B is the portion of rigid body modes corresponding to boundary points, 
iRI ^ 

ij)_ is the complementary part corresponding to internal points, 1 is an identity 
* iC 

matrix, <fr_ defines the static deformation shapes of the substructure internal to 
R IN 

its interface boundary, and represents the fixed boundary normal modes. In 
evaluating the matrix product 

1 i T i 1 
c ' ♦* " ♦* 

in this case, p* may be partitioned into submatrices corresponding to boundary 
and internal coordinates 


IBB . iBI 

p i p 

ilB 1 ill' 
P , P 


The symmetric matrix c is given 


iRR 

1 

l 

IRC 


iRN 

c 

1 

c 

1 

c 

ICR 

1 

1 

iCC 

1 

1 

iCN 

c 

1 

c 

1 

c 

"iNR“ 

f 

' Inc 

T 

"inn" 

c 

1 

c 

i 

c 


where, temporarily dropping the superscript i and the subscript R for notational 
convenience, 


RR ,RB T BB .RB . ,RB T BI ,RI . ,RI T IB ,RB . A RI T II A RI 

p <|> + 4> p + $ p <j> + $ p <f> 

RC ,RB T BB .RB 1 BI ,C ,RI T IB . ,RI T II ,C 

c * p + <f> p $ + $ p+$ p $ 

RN RB T BI .N . RI T II ,N 

c * <J> p 4> 4- <j> p <f> 

CC BB . BI A C . X C T IB . a C T II X C 

C *p + p <j> 4* $ P + ^ p 4> 


(13a) 

(13b) 

(13c) 


CN BI ,N , .C II .N 

c ■ p $ + f p $ 


NN .N II .N 
c * <|> p <)i 



I 
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NN 


Only c given by (13f) may be obtained directly from fixed boundary vibration 

RR RC RN CC CN 

testa. In general, c , c , c , c , and c will not be completely determined. 

There are various ways of assigning values to these submatrices. It is 
first recognized that a representation of p**may be derived by the following 
method. By definition, 


Given that 


it follows that 


N II . N NN 
!> P * * c . 


. N II .N 
(> U ♦ I, 


II 


, N J J 

In particular, (14) will be satisfied if 


N II . 
♦ P * 


N | NN( N T II N ) 
/c \4f u * / • 


(14) 


II 


.II 

P 

II 


II ,N NN , N II 
v $ c v • 


(15) 


It is to be emphasized that p is not determined uniquely by this method. In 
fact, may not even resemble p**. Nevertheless, 0** 
will satisfy 


NN 


»* e 11 ♦* 


and one may determine values for the submatrices (13a) through (13e) by neglecting 

all but the last term in each equation and substituting p** « p**. This is equiv- 

BB BX 

alent to assuming that p and p are both null. If this assumption it not 
acceptable, then additional tests will presumably have to be made. 

This line of reasoning was not pursued much further because of practical 
difficulties associated with obtaining the boundary modes themselves. These are 
discussed more fully in Reference [7]. From the standpoint of damping, however, 
it is apparent it the outset that both the free-interface mode method and the 
fixed-interface mode method suffer from the same basic limitation: the repre- 

sentation of damping near interfacial boundaries will tend to be poor. In the 
free-interface mode case, the tendency is caused by the fact that convergence of 
system modes is found to be poor, with more substructure modes required to attain 
a given degree of accurac". This is usually explained by the failure of lower 


frequency modes to edequately "work" the local structure near Interface boundaries. 
If strain energy in this region is not properly accounted for, it is likely that 
the dissipative energy may not be either. 

In some respects, the problem associated with using fixed interface 
modes for damping synthesis is similar. Displacements near interface boundaries 
are usually defined by static displacement shapes for which damping information 
is not available. So again, the basic problem is an inadequate representation 
of the structure's dissipative properties near substructural interface boundar- 
ies. 

On the other hand, the use of mass-loaded interface modes can poten- 
tially overcome this problem. Mass loading causes the structure to be "worked" 
more in these local areas. It is felt that the dissipative properties of the 
structure can therefore by represented more fully. Selection of the proper mass 
loading will be important. This topic is also discussed in (71. Since the 
damping matrices in this case are of the same form as in the free-interface mode 
case, their evaluation is straightforward. 
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6. Use of 


>ed Modes to Compute the Modal 


>ing Matrix 


In searching for a way to more fully describe the damping properties 
of a structure, two things become apparent. The first Is that most structures 
tend to be lightly damped, therefore damped modes closely resemble the hypothe- 
tical undamped modes except for small differences In phasing. The second Is 
that phase separation techniques presently employed In vibration testing are 
potentially capable of yielding these small ohase differences qppmtitetiyel^. 
This phase information can be used to derive a more complete set of structural 
damping Information than simple decay tests provide. 

Intuitively, one would expect that If the damping in a structure could 
be gradually reduced to zero, the damped modes wr ' 1 approach the undamped modes 
In some continuous fashion. Thus, a damped mode • ^ht be representable by a lin- 
ear perturbation of the corresponding undamped mode. This thinking led to the 
formulation of the perturbation anal,, ris contained in Reference [6]. The salient 
features of the derivation are included here using somewhat different notation to 
pro.-ide continuity in the overall formulation. 

The basic equations of motion for a substructure are first considered 
at the discrete coordinate level, as given by (1). Transformation from the dis- 


crete coordinates x 1 to distributed coordinates p 1 using x* 


♦ R p* results in 


i~i , i*i , .1 
mp + c p + k p 


where 


T 

i ,i A i i 

® - <t> R u x 4> R 

i i T li 

C " *R p *R 

k 1 - 

k - * R *c * R 


(17a) 


(17b) 


(17c) 


£ p " *R f (17d) 

The modal damping matrix c* may assume different forms, as suggested by (10) and 

iNN 

(12). The main objective is to determine the submatrix c which appears in 
both equations. For convenience, (10) and (12) may be put in the same form by 
writing both (9) and (11) as 


. i [\iB| IN 
R ■ |_*R , ♦* _ 


J 



J 


J 


1 
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With this sort of partitioning , (16) becomes 


m lBB ‘m iBN 

to iB ) 

IBb! iBN 

r I p 


. iBB ' iBN 
k 'k 

W 5 

( iB) 

In ff 

f iB 

iN 

"iNB 7 Inn 

m ,m 


v C 

”inb t Inn 

c ,c 


) In” t " | 

(p ) 

- J 


m «■ 


. - 

| 

P 1 


Then one may write 


where, in general. 


INN” IN . INN* IN , .INN IN ;iN 
m p + c p +k p ■ f 


(18) 


(19a) 


s iN ,iN 1NB-1B iNB'iB , iNB iB 

f - f - m p - c p -k p 

P P 


(19b) 


Now it is recognized that when 


A iN i A iN _ _ 

♦r “ ♦* = 1 

IN 

by proper normalization of , that (19a) is of the form 


( 20 ) 


Iq + 5q + q - f q 


( 21 ) 


where u>^ Is a diagonal matrix of undamped frequencies squared, and £ may be 
viewed as the modal damping matrix. 

From this point on, it will simplify the notation considerably to drop 
the superscript notation used to distinguish different substructures and differ- 
ent classes of modes for a given substructure. It will follow, for example, that 

1 i . ^iN 

x - x, \i - m, and <J> R * $ R < 

form of Equation (1) becomes 


x* * x, m* * p, and With the simplified notation, the homogeneous 

K R 


•• « 


mx + px + tcx ■ 0 


(22) 


These second order equations may be written in first order form leading to the 
eigenproblem 


([**] • 


'J 


M • 0 
I 


L .) 


(23) 


where A^ is a complex eigenvalue of the form 


W ‘“3 


(24) 



1 


J 


] 


_ . . J 
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© 


and $ is that part of the complex eigenvector which corresponds to the displace- 
ment vector x. The term -o. ■ c.u may be interpreted as the decay rate associ- 

3 3 J 


ated with mode J, while w 


J 


U) 


; , c 2 i® the damped natural frequency of that 
°j 

mode. The vector will henceforth be referred to as the jth damped mode . It 
may be related to the hypothetical undamped mode by the equation 

R j 


+ 5$ + i6q> 

1 Kj R j 


(25) 


where 6$ and 6$ T are considered to be small compared to . For lightly 

*j T J R 1 

damped structures governed by Equation (22) , the eigenvectors will occur in con- 
jugate pairs. The resulting matrix of eigenvectors will transform (23) to diag- 
onal form (assuming that u, p, and k are symmetric) so that, in particular, 


T 1 T 

j Jt. 
$* T 1 x*** T 


I 

_p_ i_u_ 

W i 0 


-♦-Lt*. . \M°- 

<|>X , <|>*X* 0 , A* 


(26) 


where the asterisk denotes the complex conjugate. Three separate equations may 
be derived from (26) 


♦X + ♦! «* ♦! < + v] ^ < 


all j, k 


+ Vk + ; j*i j i 


(27a) 

(27b) 

(27c) 


Using Equations (24), (25), and (20), and the perturbation assumptions 

| 0j |< < Uj 

' l* Rj l 


I pf 


r j I * < I 'TV* ^ I 


O 


one may derive from (27a) the relationships 

.T 


MJ 


♦ij 0 ♦tj ' - 2o j - 2c 3 “oj 


<jk ■ \ " \ " - ( » J -“k )( +R j “ s h k - { *i j ; J * k 


(28a) 


(28b) 



] 


1 


1 
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From (28b) and (27b), an alternative expreaaion for £,. may be found. 


T T 

« 4 <54 T u 4„ + «!,♦« W 6$, 


"Jk w J^I j - T R k 
Finally, from (27c), one may derive 


- 2 *^- 2 ^ +1 j 


a relationship which will be needed later. 

It will be of interest to examine more closely the structure of 
Equation (29), and compare it to that of Equation (28a). It has been postulated 

that £.. and £ are of the same order of magnitude (at least when u > to ). 

J* JJ °i °k 

This Implies that the elements of 64- should be of the same order J 

A A 

as 6. . Since ] 6^ )< <|$ ) ,, the J phase angles associated with elements 

J 1 

of the complex vector 4j J will be of the order 6^. This relationship 
will govern the accuracy requirements on phase angle data. 

The Important conclusions to be drawn in this section are that the 
diagonal elements of the modal damping matrix £ correspond to the uncoupled 
modal damping constants for each mode, and that the off-diagonal elements can 
be determined provided that some way can be found to identify 6$^ . 


I 


1 


I 


! 

1 
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7. Use of Resonant Response Data to Identify Dapped Modes 

The method for constructing a nondiagonal damping matrix outlined in 
Section 6 requires that the damped modes of the structure be known. These modes 
are complex, having both real and imaginary parts, and occur in conjugate pairs 
for lightly-damped structures. In a practical sense, it is not likely that these 
modes will be directly measurable from vibration tests because it is difficult to 
excite sufficiently "pure" modes. Off-resonant modes may contribute signifi- 
cantly to the total accele. ^cion response, particularly in th< ‘incident com- 
ponent defined to be in phase with the forcing function. The .rature compon- 
ent leads the force by 90 s and is normally used to define the leal undamped 

modes. The fact that off-resonant mode rerponse tends to be relatively small 
and out-of-phase with that of the resonant mode has provided the basis for phase 
separation techniques now used to more accurately define the undamped modes; how- 
ever, contamination of the coincident response by off-resonant mode participation 
imposes the need for modal separation if the damped modes are to be determined at 
all. An Iterative procedure has been developed and demonstrated. For damping 
levels on the order of IX, convergence has been achieved in only one or two 
iterations although the process may be cycled any number of times should con- 
vergence proceed more slowly. 

The basic technique to be used for mode separation requires, in addition 
to resonant response frequencies and damping rates, a knowledge of the force 
input to the structure. Assuming that the quadrature component of each reson- 
ant response is a reasonably accurate representation of the undamped mode associ- 
ated with that frequency, the extent of off-resonant mode participation can be 
determined to good approximation and subtracted from the total response, leaving 
an improved representation of the resonant mode. This procedure was successfully 
applied by Stahle [8] a decade ago to Improve upon the real undamped modes of a 
structure. The present work extends that method for application to the damped 
modes which are complex. The generalization is valid provided that the struc- 
ture is lightly damped. Specific validity criteria which involve quantitative 
relationships among the distribution of the forcing function, the amount of 
damping, frequency separation of the modes, and measurement accuracy are derived 
in the Appendix. In keeping with the simplified notation, the equations which 
follow may be considered to apply to a substructure even through the superscripts 
are dropped for notational convenience. | 

i 

i 
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It le assumed that the n equations of motion for a real structure 
(substructure) may be written In the form 


pX -f px + kx ■ f(t) (31) 

In order to make this development applicable to free-free ntructu ee as well as 
constrained structures. It will be assumed that the displacement vector x defines 
elastic deformations relative to rigid body motion. In the free-free case, this 
Implies that to utilize total response measurements, rigid body response must 
first be computed independently and subtracted from the total response. This is 
discussed further In Section 9.5 

Under these conditions, the equations of mrtion as given by (31) may 
be transformed accoidlng to 

x - 4 » r q (32) 

where is a modal matrix containing undamped elastic modes. Then (31) becomes 

Iq + 5q + u>* q - f(t) - f q (t) (21) 

The force vector f(t) will include any forces which couple the elastic modes 
to other modes, as indicated by (19b). This equation was given earlier as in- 
dicated. it ie assumed that, in general, a different force vector f(t) will be 
generated to excite a resonant response at each natural frequency , and conse- 
quently, that 

f(t)-P x 8j (t) (33) 

The complex frequency response vector corresponding to is then given by 

■> 


' [(“= - B j \ l ) + v] \ 


Initially, it may be assumed (as shown in the Appendix) that 


s jk “°Jk^ c 


+ °k) 


where 6^ is the Kronecker delta, and that is equal to the normalized quadra- 
ture component of acceleration response. Thus, a vector Hq^ may be evaluated for 
each of the resonant response condition? and transformed back to the x coordinate 
system by 
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Combining these equations in a single matrix equation gives 


I I 


H 


♦A 


where is assumed to be a nonsingular square matrix. Then 


♦« 


H H -1 
x q 


(35) 


provides the first estimate of the undamped modes . The vector H X j is obtained 

from total acceleration response data as shown in the Appendix. In this way, an 

improved version of $ is obtained. 

R 

In order to derive an estimate of 6$^., Equation (31) is written in the 
first c .der form 


_P_ l _M_ j_ x _| + k ^ |_ x 

y i 0 | x ) Oj-ufx 


.m _ 
0 


Solution of the corresponding first order eigen-problem leads to the trans- 
formation 

( 


x 

X 


p- H 

. 4 X J ( ) 


♦x ; * 

The frequency response vector in the z coordinate system is then given by 


Hz j “ 


(in i-xy 

- - J - --k - - - - 

”1 

r* i 
- - a- 

0 

-1 

♦ T 

o '(in, i-x*) 
i j 


0 1 
1 

A* 


<|,* T 


x j 


(36) 


where A is a diagonal matrix whose elements are given by (30). If it may be 
initially assumed that 

Aj * i2u)j 

and that <j> is given by the real part ot the resonant response matrix, then a vec- 
tor H 2 j may be evaluated for each resonant frequency and transformed back to the 
x coordinate system by 

H - [ ♦ • ♦*] H 

j j 

Combining these equations into a single matrix equation gives 


H 


t ♦ i ♦*] v 


(37a) 


i , 
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Alternatively, recognizing that $ and may be expressed in terms of 
$ R + 5 *R * and one “a? write 


“* ' *R \ + 5 *I H I 


(37b) 


where and Hj are both complex square matrices. Having already found $ R from 
(35), 6$^ is given by the equation 

5 *I ’ - ♦* V "I 1 (38) 

With $ and 6«p T computed, the first iteration is complt '-*. To begin 
the second iteration, these values are used in (29) and (30) to update the 
matrices £ and A, respectively. Then Equations (34) through (38) may be recycled 
to complete the second iteration. Succeeding iterations follow in the same way. 
The entire procedure may be summarized by the following: 

Initial values 

Define the damped modes initially by the quadrature acceleration response 
QUAD (h»). That is, 

(0) = QUAD (H.. ) / [QUAD (H-. ) T p QUAD (H„ )] 1/2 

j X j X j X j 

These vectors are defined to be real, and are of the form <J> ■ <J> ;• the imaginary 

J Rj 

parts 6$ = 0. 

J 

Evaluate matrices C and A 


(1) 5 )k ‘ _S jk ( °j + V ‘ <“j ' V (*Rj * 


f T T 

i x k 1 


j ’ \) 


(2) A. - Vi + i 2 oj 

J 3 R j j 3 


where 6.. is the Kronecker delta 
jk 


. i 1 • J 

l 0 : j 


J 
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(6) H x - ( * ; ♦»] H s = * R \ + «» R Hj 

(7) «» t - (h, - * R Hg) H^ 1 


Steps (1) through (7) may be repeated any number of times. The procedure should 
converge if the proper conditions are met. These conditions are derived 
in the Appendix. 


Note that P is used in Step (3), while P (without the bar) is used in 
Step (o). ^ 
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8. Analytical Verification 

A variety of examples have been worked in the process of formulating 
equations and coding them for numerical computation. They were used as check 
cases to debug computer codes and to facilitate the understanding of unfamiliar 
concepts. As such, they have been of considerable value in relating the theory 
to practical application. Some of the examples are included in this section. 

The first set of examples lv be discussed in Section 8.1 was selected 
by GAC to demonstrate the relative merits of coupling procedures involving sub- 
structure modes with different kinds of boundary conditions. The purpose of that 
investigation was to provide a basis for selecting a coupling procedure best 
suited for application to the Shuttle. Since the matrix method for damping synthe- 
sis relies on the same transformations as the synthesis of undamped equations, it 
was a simple matter for GAC to punch them on cards and transmit them to TRW. This 
avoided duplication of effort and provided a consistent set of examples by which 
to evaluate the convergence of modes and frequencies as well as modal damping in 
representative Shuttle configurations. 

In Section 8.2, a four-degree of freedom lumped parameter model is intro- 
duced to demonstrate the computation of a modal damping matrix. In this example, 
the damped modes are first computed by a complex eigenvalue routine. The modal 
damping matrix is then computed without using the original damping matrix, using 
only the mass matrix, the frequencies, and the damped modes. For comparison, the 
standard computation is made using the known damping matrix. The resulting matrix 
is nondiagonal in form, confirming that the damping is nonproportional. 

The same example is used in Section 8.3 to demonstrate how the damped 
modes can be deduced from resonant response information. In this case, it is 
assumed that neither the original damping matrix nor the damped modes are known, 
only response at each of the system’s resonances. The response includes contri- 
butions from all the modes in general, but mainly from those which are adjacent 
in the frequency spectrum. It is shown that the off-resonant mode response can 
be removed, leaving just the damped modes themselves. Computations involve an 
iterative procedure which alternately computes the damped modes and the modal 
damping matrix until both converge. 


i 



1 


I 


.] ... - 


- «... 1 


18058-600 1-RU-00 
Page 25 




8. 1 A Numerical Comparlaon of Coupling Procedures 

A comparative evaluation of the coupling procedures referred to in 
Section 5 has been made on the basis of math model data furnished by GAC [7]. 

These data reflect the mass and stiffness properties of realistic Shuttle con- 
figurations. Damping parameters for these models have been assumed arbitrarily. 

The main objective has been to define simple nonproportional damping matrices 
for each substructure corresponding to its respective discrete coordinate system. 

Modal damping matrices are then derived from the substructure mode matrices pro- 
vided by GAC. These are treated in the same way that experimentally derived 
matrices will be treated in subsequent work. 

The Shuttle configuration adopted in this study involves two substruc- 
tures, orbiter and tank, which are connected in a parallel arrangement. Schematic 
drawings representing each substructure are shown in Figure 1. The coupled con- 
figuration is shown in Figure 2. The models are planar, having three coordinates 
defined for each nodal point: two translations, and one rotation. No mass is 
associated with any rotational coordinate. Thus, a total of 14 component modes 
are obtained for the orbiter and 22 for the tank. Pinned connections are assumed 
at the attachment points. With six constraint equations acting to couple the two 
vehicles, the total number of degrees of freedom for the coupled vehicle totals 30. 
Eigenproblem solutions including all 30 degrees of freedom are used for reference. 
Truncated mode solutions involving 16 and ten degrees of freedom are obtained to 
evaluate convergence. The former is understood to include 11 uncoupled component 
modes for each substructure, while the latter includes eight, the lowest in each case. 

Numerical computations for four different coupling configurations have 
been made. These configurations are: 

Case 1. A free-free orbiter coupled to a free-free tank. 

Case 2. An orbiter fixed at its junction points coupled to a 
free-free tank. 

Case 3. A free-free orbiter coupled to tank fixed at its base 
whose junction points are free. 

Case 4. An orbiter fixed at its junction points coupled to a 
tank fixed at its base whose junction points are free. 
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Cases 1 and 2 each include two subcases designated by a and b. Cases la and 2a 
involve the use of free-interface substructure modes when the substructure is 
considered to be free-free. Cases lb and 2b correspond to substructure modes 
reflecting mass loading at the junction points. 

Computations are executed in the following manner. 


1. Preset Cjj ■ 0, all j. 

i fT j .» 

2. Read T and evalute + Tj c T* for all j. 

3. Read p* and evaluate c* * <J>* p* 4>*. 


4. 


Read $ 


i 

R * 


The procedure is recycled until all (both) of the substructures are taken into 

account. Matrices <(>* and T were furnished by GAC. 

R i 

Several different sets of p damping matrices were assumed, the simplest 
corresponding to externally-grounded lateral dashpots, as shown in Figure 3. 
Although this damping distribution gave diagonal P* matrices, the corresponding c* 
matrices were not diagonal. In order to specify internal damping, the damping 
matrices must involve rotational as well as translational coordinates since coupled 
modes may involve rigid body rotations of the two substructures. Such damping 
matrices may be derived in the same way as the stiffness matrices. Although this 
form of damping is more representative, it is not deemed essential for evaluation 
of coupling procedures, and time constraints were prohibitive. Therefore, only 
external damping was assumed. This permitted all rotational displacements to be 
deleted from the substructure modal matrices 4>*, thereby reducing somewhat the 
computational effort. 

The nominally chosen distribution of damping is illustrated by Figure 3 
where dashpots are shown connecting each substructure to ground. The dashpot con- 
stants were assigned identical values of 300 lb-sec/ in. Tabulated results from 
the damping synthesis involving Cases one through four are given in Tables one 
through four, respectively. Damping constants are given for seven elastic 
modes in each case. 
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Tigure 3. Nonproportional Damping Model for Orb iter and Tank 










b. Mass-Loaded Junction Points in Free-Free Modes 


1 

4.295 

E-3 

6.565 

E-3 

6.840 

E-3 

2 

9.309 

E-3 

7.831 

E-3 

7.644 

E-3 

3 

1.335 

E-3 

2.627 

E-3 

2.630 

E-3 

4 

1.578 

E-4 

1.530 

E-3 

1.540 

E-3 

5 

5.942 

E-5 

1.073 

E-6 

2.279 

E-6 

6 

2.851 

E-6 

6.628 

E-5 

4.029 

E-4 

7 

1.878 

E-7 

1.169 

E-3 

6.095 

E-4 


Table 2 - 


Damping Constants C By Modal Synthesis for Case 2 
a. No Mass Loading in Free-Free Modes 


1 

7.494 

E-3 

9.424 

e-a 

6.666 

E-3 

2 

2.692 

E-3 

9.753 

E-3 

7.822 

E-3 

3 

2.054 

E-4 

1.264 

E-3 

2.630 

E-3 

4 

1.267 

E-3 

1.187 

E-3 

1.539 

E-3 

5 

1.958 

E-4 

8.159 

E-7 

2.286 

E-6 

6 

7.017 

E-5 

2.176 

E-4 

4.029 

E-4 

7 

1.097 

E-5 

8.608 

E-4 

6.100 

E-4 
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b. Ma»*-Load*d Junction Point* in Free-Free Modes 


Mode No. 

10 d • o * f • 

16 d.o.f. 

30 d.o.f. 

1 

6.087 E-3 

6.843 E-3 

6.843 E-3 

2 

8.096 E-3 

7.641 E-3 

7.640 E-3 

3 

1.175 E-3 

2.632 E-3 

2.630 E-3 

4 

1.276 E-4 

1.542 E-3 

1.539 E-3 

5 

1.149 E-4 

1.302 E-6 

2.284 E-6 

6 

4.475 E-6 

2.814 E-4 

4.029 E-4 

7 

2.365 E-7 

1.167 E-3 

6.103 E-4 


Table 3 - Panging Constants By Modal Synthesis for Case 3 


1 

8.943 E-3 

6.075 E-4 

6.663 E-3 

2 

2.033 E-4 

9.526 E-3 

7.592 E-3 

3 

1.284 E-6 

7.237 E-4 

1.668 E-3 

4 

1.181 E-4 

6.847 E-4 

1.489 E-3 

5 

1.354 E-7 

5.744 E-7 

1.890 E-6 

6 

2.052 E-7 

3.135 E-5 

3.955 E-4 

7 

3.188 E-10 

8.699 E-5 

1.618 E-4 


Table 4 - 

Damping Constants C.. 

By Modal Synthesis 

for Case 4 

1 

3j‘ 

9.321 E-3 

9.035 E-5 

6.843 E-3 

2 

3.347 E-4 

9.613 E-3 

7.640 E-3 

3 

5.630 E-6 

7.560 E 4 

2.630 E-3 

4 

1.141 E-3 

1.107 E-3 

1.539 E-3 

5 

1.378 E-6 

7.544 E-7 

2.283 E-6 

6 

2.198 E-6 

2.127 E-4 

4.029 E-4 

7 

4.490 E-ll 

8.973 E-4 

6.033 E-4 
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It nay be obaerved from these results that all but Cases lb and 2b 

yield very poor results for all modes. In Case lb, none of the modal damping 

constants have converged in the 10-d.o.f. approximation, while four of the 13 

elastic mode damping constants appear to have converged reasonably well in the 

16-d.o.f. approximation. Case 2b, based on fixed-interface orbiter modes, is 

the better of the two. Reasonably good approximations are indicated for the 

first two elastic modes in the 10-d.o.J. case, while excellent approximations 

are given for the first four elastic modes in the 16-d.o.f. approximation. 

Although the use of fixed-interface substructure modes produces the 

best results from the standpoint of modal convergence, damping associated with 

BB 

the boundary modes is difficult to assess. In previous examples, both P and 

Bi 

P have been null. For comparison , Case 2b was evaluated with lateral dampers 

BI 

connecting the orbiter attachment points to ground. Sul matrix P remains null 
BB 

in this case, but P does not. Damping constants C^ for the first four elastic 

modes of the 30-d.o.f. solution are compared in Table 5 to those obtained prev- 

iously for P * 0. Considering P i 0 to provide the reference solution, one 

BB BI 

observes that having to neglect P (and P ) destroys some of the advantage 
offered by the use of fixed-interface modes as far as damping calculations are 
concerned. Intuitively, this effect should diminish, however, as the ratio of 
internal coordinates to boundary coordinates increases. 

In all of the^e examples so far, full modal damping matrices c* have 
been utilized. To show what happens when off-diagonal terms of c* are neglected. 
Case lb was evaluated on this basis. Results corresponding to the 30-d.o.f. sol'- 
tion are presented in Table 6 for comparison to the previous case. In this case, 
all of the off-diagonal terms were neglected, not just those of c iN *l Had only 
those of c been neglected, the comparison illustrated in Table 6 might have 
been closer. 

The damping in Mode No. 5 is quite small, presumably because this is 
a longitudinal mode whereas the others are primarily lateral and no longitudinal 
dashpots have been included up to this point. The addition of one longitudinal 
dashpot (300 lb sec/in) to the forward ends of both Orbiter and Tank produced the 
results shown in Table 7. Although damping in the 5th mode is shown to increase 
appreciably on a percentage basis, it remains small compared to the other modes. 
This may be explained by the fact that there are only two longitudinal dashpots 
compared to ten in the lateral direction, and the likelihood that the modal mass 
is larger in this mode due to essentially rigid body motion of the Orbiter and 
Tank in opposing directions. 
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BB 

Table 5 - Effect of Neglecting p In Modal Panning Computations 

Case 2b 


Mode No. 


.BB 


Damping Constant ,_C^ ^ 


p b No 


1 6.843 

2 7.640 

3 2.630 

4 1.539 


E-3 

7.165 

E-3 

E-3 

7.857 

E-3 

E-3 

2.707 

E-3 

E-3 

2.050 

E-3 


Table 6 - Effect of Neglecting Off-Diagon al Terms of the Modal Damping Matrix 

in Synthesis of o 44 
755~e IB 1J 

Damping Constant. _ 



Off-Diagonal Terms Neglected 

Reference 

1 

5.782 E-3 

6.840 E-3 

2 

3.930 E-3 

7.644 E- j 

3 

1.577 E-3 

2.630 E-3 

4 

1. £.95 E-3 

1.540 E-3 

5 

4.755 E-6 

2.279 E-6 

6 

2.006 E-4 

4.029 E-4 

7 

3.834 E-4 

6.095 E-4 


Table 7 - Effect on C. . of Adding Longitudinal Dampers to O rbiter _and Tauk 30-d.o.f. 

Case lb 

Damping Constant. 

Longitudinal Dampers Added Reference 


1 

8.352 

E-3 

6.840 

E-3 

2 

8.599 

E-3 

7.644 

D-3 

3 

2.707 

E-3 

2.630 

E-3 

4 

2.053 

E-3 

1.540 

E-3 

5 

8.278 

E-6 

2.279 

E-6 

6 

5.416 

E-4 

4.029 

E-4 

7 

7.622 

E-4 

6.095 

E-4 
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8.2 Approximation of the Modal Damping Matrix 

This section is included to demonstate the validity of the perturbation 
results presented in Section 6. The modal damping matrix of the spring-mass sys- 
tem shown in Figure 4 is sought. In this case, the free vibration of the system 
is of concern so that the input g(t) shown in the figure may be Ignored. The 
equations of motion are of the form 


where 


px + px + KX * 0 


.5 


.04 -.02 
.03 

Sym 


0 

,01 

.02 


0 

0 

-.01 

,01 


Sym 


0 

-1 

2 


0 

0 

-1 

1 


The eigenvalue problem given by (23) was solved to obtain X and $ as defined by 
(24) and (25), respectively. Numerical values are given in Table 8. The modal 
damping matrix c (which is identical to £ in this case) was approximated using 
(28a) and (29), 


• o 


c jj " c jj = ~ 2a i 

■ Sk ■ “j “ \ + \ v s \ ; 3 * k 


(28a) 

(29) 
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J 
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and also computed directly by (17b), 


♦» 


(17b) 


assuming that $ « $ + 6$ to first-order approximation. The results are com- 

R R R 

pared in Table 9 where they are shown to be in good agreement. 


Table 8 - Eigenvalues and Eigenvectors for a 4-d.o.f. System 

(a) Eigenvalues 


Coordinate 

Number 


J ; 


Mode No. 
1 
2 

3 

4 


Real (q . ) 
-2.083 E-3 
-1.062 E-2 
-2.579 E-2 
-3.317 E-2 


4.681 E-l 
1.236 
1.745 
2.557 


(b) Eigenvectors 


Real (* R + 6* R 


Mode 

1 

Mode 2 

Mode 

3 

Mode 

4 

3.503 

E-l 

-6.833 E-l 

6.382 

E-l 

-5.430 

E-2 

6.238 

E-l 

-3.233 E-l 

-6.675 

E-l 

2.467 

E-l 

7.606 

E-l 

5.306 E-l 

6.001 

E-2 

-1.066 


8.142 

E-l 

9.793 E-l 

6.960 

E-l 

1.108 




Imaginary (6<J> j. ) 





-4.570 

E-4 

6.809 E-3 

1.451 

E-2 

-2.476 

E-3 

-7.785 

E-4 

-4.638 E-3 

-4.778 

E-3 

2.286 

E-3 

-2.258 

E-4 

-2.167 E-3 

-7.698 

E-3 

7.952 

E-4 
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8. 3 Identification of the Damped Modes 

To demonstrate the use of resonant response data to identify damped 
modes, the four degree of freedom structure illustrated in Figure 4 is again con- 
sidered. Assuming sinusoidal excitation of the structure as shown, the accelera- 
tion frequency response was computed for each of the four masses including both 
amplitude and phase angle or, equivalently, real and imaginary components of 
acceleration relative to the forcing function. Numerical values are given in 
Tatle 10. The total acceleration amplitudes are plotted as functions of fre- 
quency in Figures 5a-d. This frequency response was then normalized so that, 
initially, the largest element of each column had the complex value (1,0). The 
real parts of the vectors were used to compute modal mass, and both real and 
imaginary parts of each vector were divided by the square root of its respective 
modal mass. The results are shown in Table 11. The normalized modal mass matrix 
is given in Table 12, where the largest off-diagonal element is noted to be 
1.37 E-2. 

Resonant frequencies were defined to be those frequencies at which the 
quadrature response of the largest displacement in a mode reached its peak value. 
Damping values were obtained from the real parts of the complex eigenvalues. 

They would be determined experimentally under practical circumstances. Given the 
force input, the Modal Separation Procedure (MODSEP) was initiated. After each iter- 
ation, the magnitude of the normalized incremental changes in both real and 
imaginary parts of the complex modes were determined from 


* k( *R > 

R j 


A k (6* ) 

i 


where the superscript is used to denote the number of iteration These results 
are listed in Table 13. For practical purposes, convergence is seen to be 
achieved in only one iteration. The improved modes, real and imaginary parts, 
are listed in Table 14. The agreement is seen to be good to one or two percent 


i 
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Table 10 - Complex Resonant Acceleration Response 


Coincident (Real Part) 


Mass No. 

Mode 1 

Mode 2 

Mode 3 

Mode 4 

1 

.1572 

-.3917 

-1.0092 

.3397 

2 

.2142 

1.2831 

- .5414 

- .7711 

3 

.0487 

1.0940 

2.9473 

1.0749 

4 

- .1578 

- .2156 

1.9115 

3.2120 


Quadrature (Imaginary Part) 


1 

32.0465 

-38.9270 

15.0036 

- 2.2961 

2 

57.0674 

-18.4483 

-15.7395 

10.5035 

3 

69.5931 

30.2452 

1.4041 

-45.5123 

4 

74.4808 

55.8524 

16.4846 

47.3419 


Table 11 - Normalized Resonant Acceleration Response 


Real Parc 

Mass No. 

Mode 1 

Mode 2 

Mode 3 

Mode 4 

1 

3.5029E-1 

-6.8291E-1 

6. 3100E-1 

-5.3202E-2 

2 

6. 2378E-1 

-3. 2374E-1 

-6.6981E-1 

2.4462E-1 

3 

7.6070E-1 

5.3055E-1 

7.4002E-2 

-1.0635 

4 

8.1413E-1 

9. 7990E-1 

7.0813E-1 

1.1134 



Imaginary Part 



1 

-2.4612E-3 

9. 5097E-3 

1.1652E-1 

-1 • 1599E-2 

2 

-3.6637E-3 

-2. 1261E-2 

-5.4412E-2 

3.4729E-2 

3 

-2. 1450E-3 

-2. 1242E-2 

-1 . 1802E-1 

-9.7437E-2 

4 

0. 

0. 

0. 

0. 


T 


T 


T 
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Table 12 - Initial 

Modal Mass Matrix 


1.0 -3.8313E-5 

4.3135E-3 

1.3142E-3 

1.0 

1. 3730E-2 

2.2449E-3 

Symmetric 

1.0 -2.8946E-4 

1.0 


Table 13 - Normalized Change in Modal Vector Length 


Iteration 1 


Mode No* 

Real Part 

Imaginary Part 

i 

3. 1019E-5 

8.1573E-1 

2 

6.8161E-4 

8.0687E-1 

3 

1.6637E-2 

9. 0720E-1 

4 

3.9757E-3 

9.8476E-1 


Iteration 2 


1 

1.2322E-5 

3. 2084E-4 

2 

6. 9856E-5 

4.0783E-4 

3 

4.9095E-4 

2.3828E-3 

4 

9.0579E-5 

3.3180E-3 


Iteration 3 


1 

2. 6201E-8 

5.4733E-7 

2 

4.2167E-7 

2 . 3002E-6 

3 

2.9780E-6 

4.7326E-6 

4 

6. 6202E-7 

9. 7453E-6 



SO H 
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Table 14 - 

■ Improved Model Vectors 



Real 

Parts 

Imaginary 

Parts 

Mode/Masa 

MODSEP 

Reference 

MODSEP 

Reference 

1 1 

3.5032E-1 

3.5027E-1 

-4. 5713E-4 

-4.5704E-4 

2 

6. 2380E-1 

6.2380E-1 

-7.7822E-4 

-7.7854E-4 

3 

7 . 6069E-1 

7.6065E-1 

-2 . 2671E-4 

-2. 2582E-4 

4 

8. 1409E-1 

8.1418E-1 

0 . 

0 . 

2 1 

-6.8342E-1 

-6.8334E-1 

6.8140E-3 

6.8086E-3 

2 

-3. 2321E-1 

-3.2328E-1 

-4. 6422E-3 

-4. 6355E-3 

3 

5. 3066E-1 

5. 3063E-1 

-2.1807E-3 

-2.1667E-3 

4 

9. 7923E-1 

9. 7934E-1 

0 . 

0 . 

3 1 

6. 3763E-1 

6. 3822E-1 

1 . 4370E-2 

1.4513E-2 

2 

•-6.6783E-1 

-6. 6749E-1 

-4. 6619E-3 

-4. 7784E-3 

3 

6. 0753E-2 

6.0011E-2 

-7.6738E-3 

-7.6975E-3 

4 

6.9664E-1 

6. 9598E-1 

0 . 

0 . 

4 1 

-5. 4558E-2 

-5. 4308E-2 

-2.4867E-3 

-2.4760E-3 

2 

2.4 .761 '1 

2.4667E-1 

2 . 2613E-3 

2. 2857E-3 

3 

-1.0657 

-1.0659 

-6.2719E-4 

-7.9520E-4 

4 

1.1079 

1.1076 

0 . 

0 . 


Table 15 - Final Modal Maas Matrix 


u <J>, 


1.0 


Symmetric 


-4 . 6318E-5 

1.0 


2. 7321E-5 
8.4953E-4 
1.0 


4. 2270E-4 
2. 2443E-4 
-4.1186E-4 
1.0 



! 


i 


i 
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except for element (4,3) which is about 20Z off. (The probable source of error 
will be discussed momentarily.) The final modal mass matrix is shown in Table 
15, The largest off-diagonal element now is 8.49 E-4, down a factor of sixty-two 
from its initial value of 1.37 E-2. The full modal damping matrix given in 
Table 16 was found to be in good agreement with the reference values except for 
fourth row elements which deviate by as much as 18.3%. This is a direct result 
of the 20% deviation in the imaginary part of the third element of mode 4. 

In the foregoing example, all four modes were used. The first three 
and first two were used in subsequent work to assess the effect of modal trunca- 
tion on £ . These results are shown in Table 17. The effect of modal truncation 
is shown to be small in this case. 

While this example problem seems to demonstrate the theoretical validity 
of the procedure formulated in Part A, it also points up a potential problem. 

The 20% error in the fourth mode, which was noted earlier, does not appear to 
be the result of a trivial input error or an error in the formulation of the 
method. Similar results were obtained for other examples; however, it was 
found that these errors could be reduced by identifying resonant frequencies 
more accurately. 

The source of the error is believed to be attributable to the fact that 
large changes in phase accompany very small changes in frequency near resonance. 
In the present example, resonant excitation frequencies were specified to five 
significant figures. Even so, changes in the fifth significant figure resulted 
in changes on the order of 100% in some of the coincident response values 
because of the small amount of damping (around 1%) . The problem is compounded 
by the fact that as damping becomes smaller, phase angles become smaller at the 


same time their derivatives with respect to frequency become larger, so that 
the sensitivity factor is of the form (d0/du>)/0. Hence, precise frequency 
determination is very important. 

In the present investigation, the resonant frequencies 0^ were assumed 
equal to the natural undamped frequencies on a one-to-one basis. In practi- 
cal applications, particularly when modes ar4 not well separated, it may be 
necessary to analytically determine modal frequencies slightly different 
from the corresponding excitation frequencies . Sweeping across would 
allow one to compute the rate of change of quadrature response with ^respect 
to frequency. The frequency at which this function is a minimum may be used to 
determine , as suggested by Kennedy and Pancu in Reference [9] . While this 
approach seeis to merit consideration, alternative procedures might also be 
sought. A need for further investigation of this problem is clearly indicated. 
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Table 16 - Comparison of Modal Damping Matrices 


Col/Row 

MODSEP 

Eo. (32) 

Eq. (33) 

1 1 

4.1663E-3 

4.1663E-3 

4.1662E-3 

2 

-1.4179E-3 

-1.4085E-3 

-1.4087E-3 

3 

-1.2800E-3 

-1.3356E-3 

-1. 3360E-3 

4 

7.4791E-4 

6.3308E-4 

6. 3325E-3 

2 2 

2. 1240E-2 

2.1240E-2 

2.1237E-2 

3 

-8.9431E-3 

-9.0619E-3 

-O.0595E-3 

4 

1.6128E-3 

1.4548E-3 

1.4544E-3 

3 3 

5.1583E-2 

5.1583E-2 

5.1581E-2 

4 

-4.2064E-3 

-4.2773E-3 

-4. 2791E-3 

4 4 

6. 6356E-2 

6.6356E-2 

6.6346E-2 


Table 17 

- Effect 

of Truncation 

on Computation of the 

Modal Damping Matrix 

Col /Row 

2 Modes 

3 Modes 

4 Modes 

1 

1 

4.1663E-3. 

4.1663E-3 

4.1663E-3 


2 

-1.4132E-3 

-1.4179E-3 

-1.4179E-3 


3 


-1. 3630E-3 

-1. 2800E-3 


4 



7.4791E-4 

2 

2 

2. 1240b- 2 

2. 1240E-2 

2. 1240E-2 


3 


-8.8476E-3 

-8.9431E-3 


4 



1.6128E-3 

3 

3 


5. 1583E-2 

5. 1583E-2 


4 



-4, 2064E-3 


6. 6356E-2 


Q 


4 
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8.4 Conclusions 

There are aeveral important conclusions to be drawn from the numerical 
examples presented in this section. The convergence of system damping values 
obtained by modal synthesis is governed by convergence of the system eigenvectors. 
Since this topic is to be discussed more fully in Reference [7], it will not be 
elaborated upon here. From the standpoint of damping synthesis, the coupling pro- 
cedures which avoid the use of fixed-bounds ry modes are preferable. The possible 
advantage of better convergence that fixed -boundary mode formulations offer is 
outweighed by the difficulty of determining both the static and dynamic modes 
experimentally, and by the difficulty of relating experimental damping properties 
to such a coordinate system. Of the nodal coupling procedures which do not rely 
on f ixed-bouniary modes, the one utilizing mass loading at Interface boundaries may 
be more desirable than the one utilizing just the free-free modes. Provided that 
appropriate means are available to determine what the mass loading should be, this 
procedure tends to converge faster with no apparent disadvantage relative to 
identification of damping prcperties. In fact, the dissipative energy character- 
istics should be brought out better for the same reasons that the strain energy 
characteristics are. 

A four-degree of freedom lumped mass model was used to Illustrate the 
accuracy attainable by the approximate method for generating modal damping matrices 
given the damped modes. For damping values on the order of 1%, it was found that 
the approximation yielded comparable results. When the damped modes were not given, 
it was shown that they coaid be derived from resonant response information by an 
iterative procedure which converged rapidly in that case. 

A significant problem was, however, indicated. As the modal separation 
procedure is presently used, it is assumed that excitation is applied to the struc- 
ture precisely at its resonant frequencies. It is unlikely that this will be rea- 
lizable in practical situations where modes are closely spaced in frequency. In 
this case, it may only be possible to excite the structure near a resonant fre- 
quency. In this case, more information will be required to identify modal frequen- 
cies to the degree of precision required by this procedure. This is an area which 
demands further investigation. 
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9. Experimental Verification 

The chief purpose of this study is to apoly the matrix method for 
damping synthesis to a real structural system and compare predicted values of sys- 
tem modal damping with those obtained by direct measurement. In th3 recent past, 

LaRC has conducted vibration tests on a l/15th-scale dynamic mcdel cf a Shuttle 
vehicle. The model consists of an Orbiter and a Booster which can be coupled by 
special spring assemblies. Each of these components is constructed of thin-walled 
tapered aluminum tubing. Propellant tanks are mass simulated by strapping lead 
ballast to each component at appropriate locations, 

Two basic modifications were made for this investigation. Previously, 
the coupling configuration included only two spring assemblies, making the inter- 
face forces between the Orbiter and Booster statically determinant. For these 
tests, a third spring assembly was added between the original two, creating redundant 
load paths. Also, pin connections were added to the spring assemblies, eliminating 
the requirement for slope continuity at the connection points. 

Since earlier tests were concerned mainly with the undamped charact- 
eristics of the structure, the inherently small amount of damping in the system 
(less than 0.5%) did not pose a problem. In this study, however, it was esti- 
mated that modal damping on the order of several percent would be required in 
order to extract meaningful data witn the available analog test and data reduc- 
tion equipment. Thus it was decided to attach three externally grounded dampers 
to the Orbiter. For comparative purposes, no damping was added to the Booster. 

It was felt that this lopsided distribution of damping would represent a worst 
case in some respects and help to elucidate significant trends in the data. 

9.1 Description of Tests 

Schematic drawings of the l/15th-scale model Orbiter and Booster sub- 
structures are shown in Figure 6. The coordinates defined in this figure corre- 
spond to the accelerometer locations used during vibration tests. The spring 
assemblies were included as part of the Booster substructure. Interface points 
between the Orbiter and Booster are defined to coincide with pinned connections 
between the two. A lumped mass distribution was defined for the system taking 
into account the physica 1 distribution of mass and requirements for accelerometer 
placement. Distributions for both the Orbiter and the Booster are given in 
Table 1J. 
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Table 18 - Mass 

Distribution for the l/15-3cale Dynamic Model 

Coordinate 

Number 

of a Shuttle Vehicle 
ORBITER 

Station x Station v 

Weight 

lbs 

1 

-0.50 


0 

6.40 

2 

18.6 


0 

62.85 

3 

28.0 


0 

4,80 

4 

39.6 


0 

1.53 

5 

46.5 


0 

9.34 

6 

65.2 


0 

3.87 

7 

77.75 


0 

1.33 

8 

28.0 


0 

74.55 

9 

46.5 


0 

10.38 

10 

65.2 


0 

5.20 

1 

0.0 

BOOSTER 

0 

1.20 

2 

15.0 


0 

83.04 

3 

20.3 


0 

11.35 

4 

32.0 


0 

48.25 

5 

38.7 


0 

6.87 

6 

45.3 


0 

48.75 

7 

57.5 


0 

8.83 

8 

66.0 


0 

21.4 

9 

105.6 


0 

32.1 

10 

134.9 


0 

17.89 

11 

20.3 


9.1 

2.2 (15.32)* 

12 

38.7 


9.1 

2.2 (14.08)* 

13 

57.5 


9.1 

2.2 (16.64)* 

14 

20.3 


0 

95.59 

15 

38.7 


0 

103.87 

16 

57.5 


0 

30.23 

17 

134.9 


0 

49.99 

18 

20.3 


9.1 

2.5 (15.62)* 

19 

38.7 


9.1 

2.5 (14.38)* 

20 

57.5 


9.1 

2.5 (16.94)* 


*with mass loading 
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Discrete dampers were connected to the Orbiter at the locations 
indicated in Figure 6. These dampers were identical, each consisting of a wooden 
ball attached to the Orbiter and placed in a glass jar fixed to ground. The 
jars were filled with Dow Corning damping fluid having a viscosity of 20K centi- 
stokes. The balls were immersed in the fluid. Relative motion between a ball 
and jar was permitted in both vertical (in and out of the jar) and horizontal 
directions. The damping characteristics of a typical ball damper were evaluated 
by rigidly attaching the ball to a shaker, with force and accelerrtion transducers 
oetween them in series. The vertical and horizontal test setups are shown in 
Figures 7a and b respectively. Sinusoidal excitation was provided at peak amp- 
litudes ranging from 0.1 g to 2 g, over selected frequencies from 25 to 500 Hz. 
Based on the assumption of linear viscous damping (damping force proportional to 
velocity), values of the proportionality "constant” c were determined from the 
peak exciting force amplitude F, the peak acceleration amplitude a, their rela- 
tive phase angle 0, and the frequency of excitation o>. 


^F sinQ 

c = 

a 


(39) 


Tabulated results received from LaRC are presented in Table 19. Damping in both 
horizontal and vertical directions was found to be linear with respect to ampli- 
tude. Only averages for the four amplitudes are shown. These values are plotted 
against frequency in Figure 8. They vary considerably over the frequency range. 
Another method for determining c is discussed in Section 9.3, where a comparison 
with these results is made. 


Substructure vibration tests were run separately for the Orbiter and the 
Booster. The Orbiter was suspended horizontally in the Orbiter /Booster plane and 
free-free modes were excited in the pitch (lateral in-plan i) direction by a shaker 
located at its tail. Each of the dampers, one at the nose, one at midspan, and 
one at the tail, are shown in Figures 9a-c, respectively. Booster viDration tests 
were made with the Booster suspended vertically. Free-free modes were excited by 
placing an exciter at the tail of the Booster. Two sets of modes were obtained, 
one without mass loading and one with. No external dampers were used on the 
Booster. 


Finally, the Orbiter and Booster were connected and suspended vertically. 
Two sets of modes were obtained by exciting the structure in the pitch direction 
at the Booster tail, then in the axial direction at the Booster tail. Three exter- 
nally grounded dampers were attached to the Orbiter as in the substructure tests. 
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Ball Da mper Test_Con Tj ^ration, VcrtlcaW fatj n. 
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Figure 7-b* Ball Damper Tes t Configuration, Horizontal Motion 
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This time, however, the Orbiter was vertical. In ord^r to introduce damping in 
the longitudinal direction and accommodate the new orientation of the Orbiter with 
respect to the Jars, different mounting brackets were used for the balls. They 
are shown for the nose, midspan, and tail of the Orbiter in Figures IGa-c. These 
brackets are stiff in both the pitch and axial directions. 


9.2 Summary of Test Results 

For each of the component and system tests defined in Section 9.1, 
the input force and acceleration response data were recorded on magnetic tape. 
The tapes were shipped to TRW for data reduction using the Spectral Dynamics 
Model SD109B Co/Quad Analyzer. Acceleration response at each station was sepa- 
rated into components in phase (Coincident) and 90° out of phase (Quadrature) 
with the forcing faction. These reduced data were recorded in digital form on 
puncned cards. A duplicate copy was given to GAC. 

In addition, measurements of modal damping were made at LaRC. Three 
methods were used based on (a) half-power-point bandwidth, (b) log decrement, and 
(c) Kennedy-Pancu [9] methods. All the modes obtained from component and system 
tests are listed in Table 20. Along with the measured frequency of each mode are 
included the damping values obtained in each case. Since most of the damping in 
the system came from the external dampers on the Orbiter, the damping properties 
of the Booster were not considered particularly significant. Therefore, modal 
damping for the mass-loaded Booster modes was measured for only a selected number 
of modes. 




I 


9 . 3 Orbiter Damp i ng Matrix 

Evaluation of the Orbiter damping matrix is considered first for 
several reasons. The Orbiter is by far the simpler of the two components. Only 
the first three bending modes were found below 500 Hz. The frequency of the 
first axial mode (although not shown) was determined to be above 500 Hz. From 
Table 20, it is noted that the frequency separation between modes is nearly an 
octave in each case. Modal damping calculations made by the three methods are 
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Hall Pamper Installation fo r Co pied Orb Iter /Booster Test 


Orb! ter Nose Damper 
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Ball Danmer Installation for Cou rU'-d Orbl te i /Booster Test 
Orblter Midspan Damper 


Figure 10-b 













Figure 10-c, Ball Damper tnstaH.it ton for Coupled Orb iter /Bop,.;* e.r Test, 
Orbltcr Tall D anin; r 
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Table 20 - Experimental M< dal Frequencies and Damping Values 



LATERAL (PITCH) 

EXCITATION OF 

ORBITER AT 

STA. 74 


Mode 

Frequency 

Force 

Percent Critical Damping* 

Description 

(Hz.) 

(lbs) 

laj 

(b) 

(c) 

1st bending 

101.40 

0.?8 

3.43 

3.33 

3.78 

2nd bending 

219.78 

0.58 

2.00 

1.90 

2.32 

3rd bending 

414.93 

0.50 

0.75 

0.78 

0.87 

LONGITUDINAL (AXIAL) EXCITATION 

OF BOOSTER AT STA. 136 

Spring axial 

145.56 

0.94 

— 

— 

— 

Spring axial 

151.51 

0.59 

— 

— 

— 

Spring axial 

162.60 

0.49 

— 

— 

— 

1st longitudinal 259.74 

0.59 

.231 

— 

.281 


LATERAL (PITCH) 

EXCITATION OF BOOSTER AT STA. 134 


1st bending 

37.93 

0.45 

— 

.290 

.281 

2nd bending 

101.62 

0.49 

.197 

.179 

.181 

spring axial 

145.56 

0.49 

.210 

.198 

.200 

spring axial 

151.05 

0.51 

.132 

.163 

.142 

spring axial 

162.86 

0.49 

.144 

.144 

.152 

3rd bending 

183.90 

0.50 

— 

(.300)** 

— 

spring pitch 

211.12 

0.50 

— 

(.200)** 

— 

spring pitch 

220.75 

0.51 

— 

(.200)** 

— 

spring pitch 

224.71 

0.51 

— 

.204 

— 

4th bending 

284.09 

0.50 

— 

.597 

— 

5th bending 

369.00 

1.00 

— 





*(a) Half-power-point bandwidth; (b) log decrement; (c) Kennedy-Pancu 
**Data unavailable. Values in parentheses were assumed for analysis 
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LATERAL EXCITATION OF ORBITER/BOOSTER ASSEMBLY AT BOOSTER STA. 134 


i 




Mode 

Frequency 

Force 

Percent 

Critical Damping 

Number 

(Hz.) 

JM 

1.01 

_ia) 

(b) 

. (c) 

1 

26.01 

2.77 

2.85 

3.01 

2 

38.73 

— 

0.96 

0.94 

0.62 

3 

57.57 

0.99 

0.76 

.49 - .96 

0.64 

4 

92.25 

1.00 

1.67 

1.37 

1.08 

5 

108.22 

0.70 

0.30 

0.256 

0.221 

6 

125.31 

1.00 

3.15 

3.36 

3.47 

7 

185.18 

0.97 

0.282 

0.366 

0.289 

8 

281.69 

1.41 

0.425 

0.412 

0.378 

LONGITUDINAL 

EXCITATON OF ORBITER/BOOSTER ASSEMBLY AT BOOSTER STA. 

135 

1 

57.43 

2.01 

— 

— 

0.563 

2 

108.3 

— 

— 

— 

— 

3 

185.52 

1.50 

— 

— 

— 

4 

256.41 

0.99 

— 

— 

0.307 

LATERAL 

1 

EXCITATION OF BOOSTER WITH MASS 
34.05 0.50 

LOADING AT 

STA. 134 
0.42 


2 

52.77 

1.00 

— 

— 

— 

3 

71.27 

0.99 

— 

— 

— 

4 

76.39 

1.00 

— 

— 

— 

r 

J 

81.49 

1.00 

— 

0.33 

— 

6 

109.76 

0.70 

— 

0.44 

— 

7 

182.81 

0.50 

— 

0.32 

— 

8 

186.56 

0.98 

— 

— 

— 

9 

187.61 

0.99 

— 

— 

— 

10 

198.41 

0.98 

— 

— 

— 

11 

202.83 

1.48 

— 

— 

— 

12 

279.32 

1.40 

— 

0.48 

— 

LONGITUDINAL EXCITATON OF 
1 53.44 

BOOSTER WITH MASS LOADING 
1.51 

AT STA. 134 


2 

65.35 

1.46 

— 

— 

— 

3 

257.06 

1.00 

— 





fjwwn t&Hr**##**-*****, **wwv' ' < * * - **f -» ; «* *j^» *V'«^«^t>i.* *+ ?* ' ^ ^ v^Tj# 


"i 

\ 


t - 1 


*r 

1 


■*•♦' '-— — * 


18Q58-6001-RU-00 
Page 64 


in very good agreement. Of course, most of the damping was placed in the 
Orbiter specifically to enhance the chances of achieving conclusive results. 

A major part of the effort was therefore devoted to avaluating the damped char- 
acteristics of the Orbiter and interpreting the results. 

This study was not without its share of problems. Certain unexpected 
difficulties did arise which caused part of the schedule to be compressed in time. 
Unfortunately, when some of the Orbiter test data appeared to be in error, it 
was not possible to repeat any tests on the damped configuration, but only on the 
undamped one. As a result, two sets of data became available, neither of which 
could be considered complete or entirely free of possible error. This situation 
required some exercise of judgement in deciding how best to utilize the data. 

There was one other complicating factor. Both sets of data led to 
poor orthogonality between the first and third modes. It was speculated that 
this might be a result of the lumped mass matrix which was used, in the sense 
that no rotational mass was assigned to either of the lead weights used to simu- 
late propellant tanks. Modal slopes are significant at these locations, particu- 
larly in the third mode. To test the hypothesis, the modes were plotted from 
available data and slopes were scaled. Additional coordinates corresponding to 
slope at Stations 18.6 and 46.5 were included, and representative mass moments of 
inertia were inserted in the mass matrix. This improved the orthogonali* y con- 
siderably. The difficulty was that modal slopes could not be scaled very accu- 
rately due to the small number of data points along the Orbiter. 

The second set of Orbiter response data corresponding to the undamped 
configuration yielded better modes than the original data since the momentum 
balance and orthogonality were both better than in the first set. The decision 
was made to use these modes in the synthesis even thongh they lacked the damping 
information. Comparison of real mode plots between the two sets indicated rea- 
sonably good agreement as far as the overall appearance of the mode shapes was 
concerned. Based on the assumption that the introduction of small amounts of 
damping should not alter the quadrature response appreciably, the quadrature 
response obtained in the undamped test was combined with the coincident response 
from the damped test to complete the second set of data. This was accomplished 


i 


i 
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by matching accelerometer response at the input locations (Station 74 in each 
case) and multiplying the coincident response from the damped tests by the ratio 
of corresponding quadrature elements at each station. The forcing functions 
assumed for the second set of data were identical to those obtained from the 
damped tests. 

In order to be as objective as possible under the circumstances, 
damping matrices for the Orbiter were computed for each of four cases: 

Case la. Original data — no rotational mass 
Case lb. Original data — rotational mass included 
Case 2a. New augmented data — no rotational mass 
Case 2b. New augmented data — rotational mass included. 

Before introducing any results, it is of interest to first examine the response 
data which comprise the basis of these computations. 

Coincident and Quadrature response at the first three resonant fre- 
quencies of the Orbiter are plotted in Figures 11 through 13. The solid lines 
correspond to the original data, while the dotted lines represent the augmented 
data. The two are in reasonably good agreement after adjusting the accelera- 
tion response at the nose of the Orbiter in the first mode of the original data 
for an apparent gain error of 10 db. It is observed that the Coincident data tend 
to display the same general shape as the Quadrature data. 

Before modal separation was begun, the rigid-body response was 
removed from the total response. Forces tending to excite rigid-body response 
included both the exciter force and the damper forces. i ^ the damper charac- 
teristics had to be known a priori. While it would have n desirable to avoid 
this by introducing some form of internal damping instead, no way of doing so was 
readily available. On the other hand, there is a distinct advantage to using 
external dampers whose characteristics are known because it provides an alterna- 
tive way of evaluating the modal damping matrix. 

Earlier when the characteristics of the ball dampers were discussed, 
an alternate method for evaluating these characteristics was mentioned. This 
method is based on Equation (28a). 
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a. Coincident and Quadrature Response 
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b. Coincident Response (Expanded Scale) 


Figure 13. Orbiter Acceleration Response at Third Resonant Frequency 
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Assuming that internal damping in the Orbiter is small compared to the external 
damping, that the three ball dampers are identical, and they they introduce 
damping forces proportional to velocity by the constant c, (28a) suggests that 



k-1,4,7 


(40) 


where the three dampers are located at coordinates 1, 4, and 7. Having measured 


values of c 


, to , and $ , values of c were computed for each of the three 

J j 


modes (Case 2b). They were not in very good agreement with those values of c 
presented in Figure 8 except for the first mode at about 100 Hz. Instead of 
c tending to decrease with increasing frequency (as in Figure 8), it increased 
as shown in Figure 14. This presented somewhat of a dilemma as to what values 
should be used. 


It was learned that 0 in (39) varied from 39-8° down to 8.9° in going 

from 100 Hz to 400 Hz in the case of horizontal motion, and from 71° down to 

14.3° in the case of vertical. The damping force on the ball was of the same 
order as the driving force and the inertial force. The "virtual mass" of the 
fluid is an uncertain factor, and could vary with frequency. On the other hand, 

the damping forces exerted on the Orbiter at resonance are much smaller than 

inertial forces, and fluid mass should have no significant effect in this case. 

It is therefore believed that values of c obtained from (40) are more reliable 
than those obtained from (39). Unfortunately, there was no clear cut way of 
checking values of c for vertical motion since, in the system tests, horizontal 
and vertical motion occur together. 

In evaluating the modal damping matrix of the Orbiter, an average 
value of c * 0.68 was assumed for the damper constant in each of the four cases. 
For comparison, a value of c s 0.36 was also used in Case 2b. 


by 


The task of computing reference values for the modal damping matrix 

8 


is not altogether straightforward either because of the variability with frequency 
of the damper "constants". However, from the form of (29), it is observed that 
in each of the two terms on the right-hand side, the frequency is associated with 




J 
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the imaginary part of the mode 5$_ , which must be related to C 4 • Other things 

T J J 

being equal, the term having the larger of the two frequencies should dominate. 

T 

In evaluating £ - p $ R , it is therefore Justified at least to some extent to 

replace the nonzero elements of the diagonal matrix p by unity, and compute a 

matrix £'. Then each element can be multiplied by the value of c corresponding 

to to > id . This sort of frequency scaling will cause the diagonal elements of 
°k °J 

the reference matrix to be identical to the diagonal elements of the other matrices, 
and will cause the off-diagonal elements to be proportioned accordingly. The refer- 
ence modal damping matrix determined in this manner is presented for comparison to 
matrices computed from the damped modes in Cases la, lb, 2a, and 2b in Table 21. 

Also included is Case 2b, computed on the basis that c * 0.36, instead of c ■ 0.68. 
This case is labeled Case 3. 

Several things are apparent in Table 21. Most noticeably, the off-diag- 
onal elements in the third column are more sensitive to differences among cases. 

This is to be expected because of the importance of rotational mass in the third 
mode. The final modal mass matrices obtained in Cases la, lb, 2a» and 2b are: 

Modal Mass Matrices 


Case la Case 2a 


"l.O 

- .097 

.623” 

fl.O 

- .010 

.543~ 


1.0 

- .123 



1.0 

- .093 

- 

Case lb 

1.0 


- 

Case 2b 

1.0 

1.0 

- .103 

.308 

1 

fl.O 

- .019 

.167” 


1.0 

- .071 



1.0 

- .043 



1.0 




1.0 




; O 


*r * i 
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Table 21 - Comparison of Orblter Modal Damping Matrices 


Case la 


Reference 


42.4 37.0 

5.77 

P 42.4 30.0 

46.8 

52.5 

16.1 

j 52.5 

18.9 


40.7 

! 

t_. 

40.7 


Case lb 

42.0 80.5 

52.5 -13.2 

40.7 

Case 2a 

34.4 35.6 

52.5 41.6 


Case 2b 


42.4 

37.2 

128.2 

f 


52.5 

8.25 




40.7 _ 



Case 3 
39.8 
52.5 


123.1 "1 


Case 2b is shown to represent the best of them. Significant improvements in 
orthogonality are seen to be offered by inclusion of rotational mass and by the 
second set of data over the first. It may also be noted that * 8 relatively 

stable and differs from the reference value by the least amount. A computation 
was also made using Case 2b, but only two of the three modes. In this case, 

£ * 39 . 0 , whicn is still in the right ballpark. 

It is of interest to look at the imaginary parts S 4 >_ of the damped modal 

x j 

vectors They are plotted in Figure 15. They do not resemble the real parts 

of the modes as closely as the Coincident response resembled the Quadrature, partly 
because of the way in which they are normalized. Ten data points were available 
to plot these curves. Eight of them are accounted for by the circled dots. The 
other two come from slopes obtained at Stations 18.6 and 46.5. Without these, it 
would have been difficult to plot 6 $ for the third mode. Each of the vectors 641 

J 1 

was normalized so as to cause tne phase angle to be zero for the largest element 
in the mode. This element happened to correspond to the tail in every case. It 
is possible that this type of normalization could obscure some otherwise meaning- 
ful characteristics of the vectors. Other forms of normalization were not 
explored, however. 


i 


i 



Damping Coefficient, c (lb sec/in) 
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9.4 Booster Damping Matrix 

Since the amounc of damping in the Booster Js very small, there was not 
much hope of computing a meaningful modal damping matrix aside from the diagonal 
terms. Nevertheless, an attempt was made to see how the modal separation program 
would work. Resonant response foe the first nine Booster frequencies with no mass 
loading was entered. The iterative procedure failed to converge. This result 
was not unexpected. In the first place, the damping is so small that the accu- 
racy of the coincident response measurements is in question, and secondly, the 
frequency separation among spring modes in particular violates the criteria for 
convergence derived in the Appendix. It was noted, however, that after just one 
iteration, the orthogonality of the modal mass matrix was improved significantly, 
in particular, with respect to the seventh mode, which is the lowest frequency 
spring mode in the pitch direction. The original seventh column of the modal mass 
matn'x is shown in Table 22 tor comparison to the same column after one iteration 
oi MOLSEP. 

Table 22 - Improvement of Orthogonality Anong Booster Mooe s 
(Column of Modal Mass Matrix) 

Row Nc ♦ Original Improved 


1 

- .116 

- .319 

2 

.055 

.134 

i 

- .094 

- .019 

4 

- .010 

.053 

s 

. *>55 

.055 


6 

- .419 

- .065 

7 

1.000 

1.000 

8 

.428 

.194 

9 

- .809 

.039 


The overall change ie seen to be for the better even though some of the elements 
did become larger. This kind of improvement was not realized for the Orbiter 
because the poor orthogonality there was caused by an incomplete description of 
„ „ the third mode, and not by having more than one mode contributing to the resonant 

■» * response. 
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The effort to compute a damping matrix was pursued further, however, by 
Introducing fewer modes in the separation algorithm. A case Involving only the 
first two bending modes was tried first. This eliminated the frequen'y separation 
problem. Convergence was achieved in a single iteration. Then the first axial 
spring mode at 145.56 Hz was Included in a three-mode case. This time, converg- 
ence (to approximately two significant figures in 6<p ^ ) was achieved in four iter- 
ations. The second axial spring mode at 151.51 Hz was added to make a four-mode 
case. Computations diverged. 

In the Appendix, a convergence criterion is derived in terms of a scaled 
impedance matrix Z(iflj) associated with each resonant excitation frequency . 

In the case of proportional damping and in the case jf diagonal mo lal damping in 
general, Z is a diagonal matrix whose largest element is unity. The other 

diagonal elements diminish ir size as they become further displaced from the unit 
element. In the case of nonproporcional damping in general, g is a fully-popula- 
ted matrix whose elements tend to diminish in magnitude a £ they become further 

displaced from the unit element in any direction. A criterion foi c mvergence is 

T~ 

that each element Z^ » e^ZCm^e^ (e^ denoting the kth column of the identity 
matrix I) of each resonant impedance matrix Z(ifij) satisfy 



where 6 ■> u> /w > 1. Some lyp^d ---..lues for Booster calculations ave given 

°k °£ 
in Table 23. 


Table 73 - Co nvergence Indicators tor Booster Damping Matri x ’- ilrulat lons 


Nu".ber of 
Modes 

hi 

. p 2 


^2; £ /(e 2 -l)‘ 


] ‘-kt ; 

2 

2,1 

7.18 

.290 

.0306 

. '583 

. 2178 

3 

3,2 

2.06 

.179 

.0581 

4. 64 

.270 

4 

4,3 

1.08 

.198 

• •• • 

* 

* 

5 

5,4 

1.165 

.li'3 

1 /,o 


* 

6 

6,5 

1.278 

.144 

. .122 

* 

* 


•Values unknown because iteration failed to converge 






The (3,2) elements of £ are seen to be quite large. This may be due to an insuf- 
ficient number of response measurements, as the results of some cf the Orbiter 
calculations suggested. It is not particularly easy to ’’el±eve that an axial 
spring mode could couple so strongly with the bending modes. What appears to 
be happening i~ that che modal mass of the spring mode is so small that the ele- 
ments of are large compared to corresponding elements of or , for 

j l 2 L 1 

example, since the *. are n rmalized to give unit medal mass. Errors would 

tend to be amplified xn this case. 

It is clear from Table 21 that the cwo— , :>de case should have converged 
rapidly since the upper bound B on was much smaller than unity. It is inter- 
esting also to note that in the three-mode case, convergence was achieved even 
though B * .270. This is not too surprising since other off-diagonal e 1 aments of 
the Z matrices are much smaller. In the cases involving four or more muucs, con- 
vergence was not achieved at all, so that values of are unknown. 

It is of interest to examine the response measurements and the first two 
damped modes computed for the Booster. Response measurements are plotted in Fig- 
ures 16 and 17, while the corresponding damped modes are plotted ir Figures 18 and 
19. Again, apparent errors were found in the original response data so that tests 
were rerun for some of the modes (including the first). Data from the second 
test are plotted with circles, while data obtained originally are plotted with tri- 
angles. Ten-db gain errors were apparent at Stations 20.3 and 45.3 in the first 
mode of the original data. 

Several things are observed in Figure 16. Coincident data from the first 
test follow the shape of the Quadrature data and are greater in magnitude than 
Coincident data from the second test, which also exhibit more scatter . Different 
procedures for data reduction were used in each case. In the first case, data 
were recorded simultaneously on tape over a 20-second time interval and reduced 




a) Quadrature Response 





Second Resonant Freouenc 






18058-6001-RU-00 

Page 81 


U using the Co/Quad Analyzer. Data from the second test were recorded manually in 

real time o.er approximately a five-minute time interval. Amplitude and phase 
information were taken, phase angles being read from a phase meter. No phase lock 
system was used in either case although preliminary tests indicated no significant 
phase shifts over a typical 20-second interval. Over a five-minute period, how- 
ever, the possibility of significant phase shifts is acknowledged. Phase drifting 
-ould explain the scatter as well as all of the positive Coincident data points. 

Another notable characteristic is that the imaginary parts of the modes 
are considerably greater in magnitude than they should be if only damping informa- 
tion were reflected (see Section 6). Yet, ” ^1^1 * A plausible expla- 

nation for this is that the mode was probably not tuned pe.iectly, causing some 
of the resonant Quadrature response to "spill* 1 over into the coincident. This 
would tend to account for the fact that the Coincident and Quadrature shape char- 
acteristics are so similar, at least in the first test. This explanation would 
also admit the possibility that computations of might not be altered appre- 
ciate^ if phase shifts were proportional to the real part of the mode. This is, 



(e 4> t 


6<t> T ) - ♦: 
h z 




1 


(where t is a small proportionality constant) because of the supposed orthogonal- 
ity of 4> and <J> . 

z R 1 

The damped mode characteristics shown in Figure 18 are interesting, too. 
Aside from the sign difference in the imaginary part of the mode, results from the 
two tests bear a fairly close resemblance. For the sake of comparison, the modal 
matrix obtained with the substitution of Test 2 data for the first resonance into 
Test 1 data is 




f 1.385 


Sym. 


.052 

2.285 


-8.382 
2.286 
3.622 . 


^ * 
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The largest change from earlier computations is seen to occur for the (1,2) 
element, supporting the contention that the highest frequency term in (29) is 
dominant . 

Although the Booster damping matrix computations are felt to contribute 
significantly to the value of this investigation, it was not deemed critical that 
the full modal damping matrix be included in synthesis calculations because of 
the small amount of damping. Since not all the off-diagonal elements could be 
obtained anyway, the diagonal matrix 


e 


1.38 


2.28 


3.62 


3.10 


2.94 


5.78 


5.31 


5.54 


5.76 


was used. No computations for the mass-loaded Booster case were made since the 
modal synthesis attempted by GAC yielded unsatisfactory results. 


* *. 

) 

+ * 


9.5 Damping Synthesis 

The modal damping matrices given in Sec Lions 9.3 and 9.4 provide only 
part of the component damping matrices c* needed to synthesize system damping. 

If the Booster is designated as Component 1 and the Orbiter as Component 2, con- 
sistent with the convention adopted by GAC, the modal damping matrix derived for 

Inn 2nn 

the Booster may be denoted by c and that for ti Orbiter by c , according to 

1 RR 1 RN 

the notation of Section 5. The other submatrices c and c must be found 
before synthesis can begin. In cases where there is no external damping, these 
matrices are null. Thus, for the Booster, 

lRR 


1RN 


0 

0 . 


2rr 2RN 

The Orbiter does have external damping, however, so that c and c 


^ IflW --rt " firhYiitrtrrii iff frrtfiWilfo 1 * 


are nonzero. 






where c and c were, respectively, the horizontal and vertical damping Constants 11 
H V 

of the ball. dampers. 

It is recognised that Orbiter coordinates 8, 9, and 10 do not correspond 
to axial displacements of the Orbiter at Stations -0.5, 39.6, and 77.75, where the 
dampers were actually attached. However, in synthesizing system modes, GAC used 
only the component modes up to about 223 Hz. Since the first Orbiter axial mode 
is above 500 Hz, axial motion of the Orbiter was primarily rigid-body motion. 

This was confirmed by the data. Small variations among these axial response mea- 
surements were present due to the positioning of accelerometers about three inches 
off the centerline of the Orbiter in the Orbiter/Booster plane. The accelerometers 
were actually mounted on the spring brackets. Orbiter bending did induce small 
amounts cf relative axial motion, but they were negligible. 


7 
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Because of the dependency of damper characteristics on frequency, it was 


hard to decide which values to choose for c H and c.^. To be consistent with earlier 
computations of the modal damping matrix, a value of c u * 0.68 was used. Figure 14 
indicates that while no corresponding values of are available from Orbiter 
damping tests, separate tests performed on a typical damper indicate that c^ >C H * 

It was therefore assumed that c^ ■ 0.900* 


The selection of these values might bo disputed on the grounds that 
according to Figure 14, c » 0.68 corresponds to a frequency of approximately 
250 Hz, while the system modes range from about 26 to 200 Hz. At 100 Hz, c^0.45. 
Maintaining the same ratio between c^ and would indicate a corresponding value 
of Cy * 0.62. Computations were made using both sets of values, and are shown in 
Table 24 for comparison to measured values. The agreement is seen to be quite 
good, especially for the case where c„ 3 0.68 and c = 0.90 

n V 



Table 24. A Comparison of Predicted and Measured System Damping 


System 

Mode 


V .45 

y - 62 


c = . 68 
h 

V 90 


c =.68* _ , 

H Reference 

Cy=* . 90 (Experimental) 



1 

2.24 

3.23 

3.21 

2.77 -3.01 

2 

.768 

1.03 

1.02 

.62 - .96 

3 

.485 

.658 

.660 

.49 - .96 

4 

.835 

.765 

.855 

1.08 -1.67 

5 

.262 

.311 

.331 

.221- .30 

6 

2.42 

2.98 

3.10 

3.15 -3.47 

7 

.193 

.197 

.219 

.282- .366 


Another set of computations using these values for c and c. , only this 

2NN ** ^ 

time a diagonal c sub-matrix was assumed for the the Orbiter by setting the 

off-diagonal element to zero. These results are also shown in Table 24, Inter- 
estingly enough, deleting the off-diagonal element made little difference. The 
reason is obvious when one looks at the modes of the coupled system. There is so 
much frequency separation between the Orbiter modes that for all practical purposes, 
they combine only one at a time in the synthesis of system modes. 


*The off-diagonal element of the Orbiter modal damping matrix was neglected 

in this case. 
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This conclusion is somewhat disappointing in the respect that a 
comparison of the predicted modal damping constants of the system to those mea- 
sured directly does not establish the Importance of including the off-diagonal 
terms. In fact, the overall effect of neglecting the off-diagonal term seemed 
to make things a little better rather than worse, especially for the fourth mode. 
One must be cautioned against drawing this conclusion too readily, however, 
because the modal damping predicted in the fourth system mode was significantly 
in error to begin with. The reason for the error is indicated by a comparison 
of the predicted and experimental modes (Reference {7]), The experimental 
mode exhibited a large displacement cf the Orbiter tail, whereas the predicted 
mode indicated practically no displacement. The predicted frequency was 98.26 Hz, 
whereas the experimental frequency was only 92.25 Hz. The system mode did contain 
a large amount of the first Orbiter bending mode whose frequency is 101.41 Hz. 


9.6 Conclusions 

Application of a matrix method for damping synthesis to real struc- 
tures has produced encouraging results. The predicted values of system modal 
damping agreed very well with measured values, for the most part. Resonant 
response was successfully used t evaluate the coupling term between the first 
two Orbiter modes. Coupling terms involving the third mode were also computed, 
but not enough response measurements were available to fully describe that mode. 
Consequently, large amounts of uncertainty were introduced in those terms. 

Damping calculations were made with the Booster test data also, even 
though damping levels were very small. Notable in this case is the improvement 
in orthogonality among the Booster modes which resulted from the modal separa- 
tion effort. While computations involving more than the first three modes 
failed to converge, the results were at least in agreement with the convergence cri- 
teria derived in the Appendix. Computations involving only two or three modes 
did converge, but there was no way of verifyirg the results. The damped modes 
contained no obvious discrepancies, and cne imaginary parts indicated reasonably 
smooth cur/es. 

The synthesis of Orbiter and Booster damping matrices to predict System 
damping yielded results wh.'ch agreed very well with direct measurement. Part of 
the success must be attributed to the use of external damping, however, so that 
the results are not as conclusive as they might be had a different model been 
chosen. 
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PART C. APPLICATION AND FUTURE DEVELOPMENT 
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10. Curient applicability 

Parts A and B document the development of a new method for damping 
synthesis based on the use of substructure damping matrices. The method depends 
on the ability to determine the damped modes of isolated substructures. These 
modes are complex, and embody phase as well as amplitude information. While the 
use of complex modes in the dynamic analysis of linear systems is not uncommon , 

It does represent a departure from the mainstream of structural dynamics analysis 
and testing. As a consequence, two major obstacles arise: conceptual difficulty 

related to the ph>sical interpretation of complex modal vectors and practical 
problems which have net previously been encountered because there has not been 
a need to acquire complex modal data. Thio section is included to help over- 
come these obstacles by providing a condensation of concepts and findings which 
are particularly relevant to practical application. 

The combined use of amplitude and phase information to evaluate dynam- 
ics characteristics is certainly not new. Classical « lalysis of feedback con- 
trol systems has frequently employed the Nyquist diagram which is a polar plot 
of output/input amplitude versus phase angle. Kennedy and Pancu [9] used the 
same representation of data for structural vibration arolysis. The primary dif- 
ference between these techniques and the present use of complex modes is that 
the former is used to completely characterize a system between d-screte points 
over a continuous frequency spectrum, whereas the latter is used to characterize 
the entire system in terms of a discrete frequency spectrum. Most current text 
books on structural dynamics describe cwo alternative methods for response anal- 
ysis: the Frequency Response Method and the Normal Mot'e Method. The iresent 
work draws upon both approaches and, in a sense, bridges the gap between them. 

The amplitude and phase response of a structure to sinusoidal excitation deter- 
mine its frequency response. This, in turn, is used to evaluate the damped 
modes. Because structural damping tends to be small, the damped modes are rela- 
ted to the classical (and hypothetical) undamped modes in a simple manner; i.e., 
to first-order approximation, the real part of a damped mode is equal to the 
undamped mode while the magnitude of the imaginary part is small compared to that 
of the real part. Through a linear perturbation of the undamped equation of 
motion, the damped modes are used to determine a full modal damping matrix for 
the structure or substructure. 
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From the perturbation equations developed in Section 6, it ia evident 
that the magnitude of modal phase angles one might expect to find would be on 
the order of c> the critical damping ratio. Thus, for a structure having 1% 
damping in a given mode, modal phase angles on the order of .01 radian, or .6 s , 
could be expected. In order to determine such phase angles experimentally, mea- 
surement errors must be relatively small. This matter is discussed more fully 
in Section 11. 

It has been pointed out that phase angles associated with the resonant 
response of a structure to sinusoidal excitation are likely to be much larger 
than C. This may be due to the presence of off-resonant mode contribution in 
the total response as well as small differences between excitation and modal 
frequencies. The effect of the latter is to "spill" some of the real part of 
the mode over onto the imaginary axis; this should be of little consequence 
however, because of modal orthogonality wi h respect to the mass matrix, as 
pointed out in Section 9.4. Off-resonant mode response must, of course, be sep- 
arated from resonant mode response in order to determine the damped modes. The 
iterative procedure developed to do this appears to work satisfactorily except 
when modes become too closely spaced, in which case, the iterative computations 
may fail to converge. Quantitative criteria for convergence are derived in 
the Appendix, and indicate that for damping on the order of 1%, modal separation 
should be at least 20% to achieve convergence. The larger the amount of damping, 
the more frequency separation is required. Damping on the order of 5% will 
require a freqeuncy separation of about 50%, in general. These conditions assume 
a ratio of unity between the respective off-diagonal and diagonal terns of the 
modal damping matrix; which is to say, the two modes in close proximity to each 
other are strongly coupled by damping. As this ratio becomes smaller, the 
requirement for frequency separation is reduced according to the square root 
of this ratio. 

Another requirement discovered in the course of this investigation is 
that thu excitation frequency at each resonance must be maintained at the reson- 
ant frequency to perhaps five significant figures when damping is on the order 
of 1%. Resonance may be defined as the frequency at which the quadrature 
response peaks for that part of the structure undergoing the largest amplitude 
of response. 
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A sufficient number of response measurements will be required to fully 
describe each mode of interest. Interpolation may be used where appropriate. 

If poor orthogonality of the "undamped" modes is evident even after modal sepa- 
ration, either a poor mass matrix or an incomplete description of the motion (not 
enough response points) or both may be suspected, assuming that the modes were 
properly tuned during test. Foor orthogonality in this case will imply errors 
in damping matrix computations since the same quantities are involved. 

Some comments are in order with regard to choosing a physical model for 
the experimental support of an effort such as this. In retrospect, it has become 
clear that the use of external damping has led to certain practical difficulties. 

In the first place, it was hard to determine the characteristics of the ball 
dampers used in this study. Different methods employed to evaluate their damping 
"constants" led to different conclusions, thus introducing considerable uncer- 
tainty. Since this information was needed to remove rigid-bc^.y response from 
the total resonant response and to account, 'or the coupling between rigid and 
flexure modes, the uncertainty had a significant influence on damping computa- 
tions. Furthermore, since the lower system modes contained large amounts of 
Orblter rigid-body motion, the influence of rlgid-bcdy damping upon predictions 
of overall system damping was strong. In this respect, the l/15th-Scale Dynam- 
ics Model of Shuttle was rot well suited to the objectives of a damping otudy. Al- 
though attempts were made to Increase the level of damping by adding internal 
damping, no method was found to provide enough damping because of the large stiff- 
ness end mass of the lead-weighted tubular structure. In the future, this problem 
can be avoided by the use of lighter weight, more flexible structures for whioh 
damping tape and the like are capable of providing sufficient damping. Flat beam 
type as well as hinged structures are possibilities for consideration. 

Of ccurse, in the case where the design of a test specimen is dictated 
by other requirements, it may be necessary to "work" with only a small amount 
of internal damping. Improvements in phase measurement accu/ary will undoubtedly 
be required. Data acquisition is discussed in the following section. 
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11. Data Acquisition 

At the outset of this study, an investigation was made to determine 
the degree of accuracy to which phase angle data could be measured with avail- 
able analog equipment [10]. Phase errors are introduced by the transducers and 
conditioning equipment, by the recording and playback equipment, and by the 
ana 1 '^ data reduction equipment. An effort was made to determine the magnitude 
of tuese errors and to minimize them whenever possible. 

Tests were conducted at LaRC to determine the variability of phase angles 
measured by diff 'rent piezoelectric accelerometers. The accelerometers were 
mounted on a rigid mass attached to a shaker. Relative phase angles between 
pairs of them were recorded versuo frequency for all 30 accelerometers used in 
this series of tests. These data fell within a + 0.1° bandwidth over a frequency 
range of 30-500 Hz. 

To assess the magnitude of phase error which could be introduced by 
recording and play-back tape systems, LaRC recorded sinusoidal signals on tape 
and shipped them to TRW for playback and analysis [11]. Static phase shifts 
tending to increase with frequency were observed (as much as 10° at 450 Hz). 

These were presumably caused by imperfect alignment of the record and reproduce 
neads. The need for phase calibration of each channel had been anticipated in 
advance, and a special calibration device was built. At tape speeds of 15 ips, 
dynamic phase variations of + 1° were observed at 450 Hz. These can be caused 
by tape stretching and low frequency flutter. It was decided to increase the 
tape speed to 60 ips to minimize flutter and operate in the tape synchronization 
mode to compensate for tape deformation. In this manner, dynamic phase varia- 
tions were reduced to + 0.1° (over a 20-second time period). 

The third possible source of error is the analog data reduction equipment 
which included a Spectral Dynamics 1012B Dual Channel Tracking Filter and a Spec- 
tral Dynamics SD109B Co/Quad Analyzer. While amplitude errors were assessed, 
phase angle error data and/or specifications were unavailable. Precautions were 
taken, however, to have the system carefully tuned for optimum performance during 
this task. 
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If the phase error Introduced by the Co/Quad Analyzer is also assumed 
to be on the order of + G.l°, the cumulative error would be about + 0.2% or 
*0035 radians. In order to get meaningful damping measurements, the amount of 
damping should therefore be greater than one percent • This was the reason for 
adding external damping to the Orbiter. 

In summary, phase errors must be held to a minimum. A phase- lock system 
is recommended to prevent the excitation frequency from drifting and intro- 
ducing phase shifts. For 1% damping, the overall phase error should be held to 
+ 0.2°. This allows for approximately + 0.1° errors in the instrumen ation, tape 
handling, and Co/Quad reduction, respectively. All data should be recorded simu- 
taneously on tape for a period of 20 to 30 seconds at a tape speed of 60 ips to 
minimize effects of low frequency flutter. Tape playback should operate in the 
synchronization mode to compensate for tape deformation. Every channel of data 
should be phase calibrated to compensate for head misalignment. 


With regard to improving measurement accuracy, a' digital system offers 
potential advantages over an analog system. By digitizing the data in real time, 
the analog tape problems are eliminated completely since records would be stored 
in digital form. The Co/Quad analysis could then be accomplished by digital fil- 
tering and Fast Fourier Transform methods, rather than having to use a tracking 
filter and the Co/Quad Analyzer. The only problem remaining would be one of 
improving the instrumentation and/or conditioning equipment. This problem has 
not been considered yet. 
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12. Comparison with Energy Method 

Two methods for damping synthesis have now been proposed. They are 
distinguished by the form in which substructure damping properties are described. 

In the matrix method they are described by a modal damping matrix, and in the 
energy method, by total dissipative energy. In the matrix method, the spatial 
distribution of damping throughout the structure is represented; in the energy 
method, a kind of average distribution is represented in the sense that total 
dissipation energy is a summation over the structure. Herein is believed to 
be the basic difference between the two methods. 

As might be expected, it is more difficult to determine the modal damping 
matrix of a structure than total dissipative energy. More data are required and 
measurement accuracy is more of a problem because of the need to measure small 
phase angles. In general, it will cost more to get the better resolution offered 
by the matrix method. 

In an attempt to compensate somewhat for the lack of resolution in the 
energy method, Kana and Huzar have proposed the use of engineering judgement to 
categorize structural modes so as to form more than one energy curve for a given 
substructure. More judgement would then be required to decide which curve to use 
when trying to establish damping energy on the basis of some given kinetic energy. 

An advantage of the matrix method is that it can be fully automated with no human 
interaction required. 

Regarding the question of linearity, it is recognized that the matrix 
method presupposes linearity in the equations of motion. This type of linearity 
was also assumed by Kana and Huzar. Although they claim this is not a fundamental 
requirement, their report does not indicate how nonlinear problems might be handled. 

One final consideration is that of uncertainly in the predicted values 
of system damping. Even when making direct measurements, it is not uncommon to 
experience errors on the order of 50Z. While the energy method is attractive from 
the standpoint that energy data are believed to be fairly reliable, it will be 
difficult if at all possible to assess errors introduced by ignoring the distribu- 
tion of damping. The matrix method is more amenable to a rigorous statistical 
evaluation of error because the information used in computations relates directly 
to measured quantities such as frequency, response and mass distribution; it does 
not involve human judgement. 
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13. Final Conclusions 

A matrix method for damping synthesis has been developed. Its 
theoretical basis relies on a linear perturbation of the undamped equation of 
motion for lightly-damped structures. Verification has been accomplished 
using both analytical models, which satisfy the inherent assumptions of linear 
viscous damping, and experimental data from tests of real structures. The 
results are encouraging in that the nonproportional damping characteristics 
of real structures have been determined for the first time, making it possible 
to synthesize the damping matrix of a structural system from those of its 
component parts. 

The present study has encompassed a number of separate investigations, 
related to damping synthesis. The following major conclusions have been 
drawn from this study: 

1. Modal damping in structural systems can be determined to first 
order approximation by operating on the substructure modal 
damping matrices with the same linear transformations used to 
couple and diagonalize the mass and stiffness matrices. 

Solution of the complex eigenproblem is not required so that 
the procedure will be easy to incorporate in existing structural 
dynamics computer programs such as NASTRAN. 

2. Coupling procedures involving fixed-interface, f ree-interface 
and macs-loaded-interface substructure modes were investigated 
for their suitability to damping synthesis. The two major 
considerations included the practical ability to determine an 
adequate description of damping at the substructure level, and 
the convergence of synthesized values of system level damping. 
Convergence of system damping depends directly on the convergence 
of system eigenvectors. While the fixed-interface mode method is 
superior from the standpoint of convergence, it requires the use 
of static modes which are difficult to determine experimentally 
and contain no damping information. The free-interface mode method 
avoids the use of static modes but converges poorly in general. 

The use of mass-loaded interface modes appears to yield a favorable 
compromise between the other two methods, offering reasonably good 
convergence without the use of static modes. Provided that a way of 
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2. (Cont’d) 

determining and applying appropriate mass distribution to the 
interface boundaries can be found, this method appears tc hold the 
most promise. 

3. The basic feasibility of determining substructure modal damping 
matrices from resonant response data has been established. The 
damping may be nonproportional, leading to a fully-populated 
modal damping matrix. Resonant response data are used to compute 
the damped modes which are also complex. These in turn are used 
to evaluate the off-diagonal terms of the modal damping matrix. 

The off-diagonal terms must in general be included in the 
synthesis. It is only when the modes of a given substructure 
participate one at a time in the system modes that the off-diagonal 
terms (while perhaps comparable in magnitude to the diagonal terms) 
may be neglected. Although this was evidently the case in the 
present application to a real structural system, it will seldom be 
true of more realistic structures. 

4. The isolation of damped structural modes from total structural 
response to sinusoidal excitation requires an iterative computational 
procedure wherein the damped modes and the modal damping matrix are 
alternately computed. Criteria for the convergence of this procedure 
have been derived and appear to be in agreement with experience to 
date. (See Section 9.4 for example.) 

5. A significant improvement in the orthogonality of the "undamped 
modes" (real parts of the damped modes) may be achieved as a by-product 
of the damping matrix computations. This will be true whenever poor 
orthogonality of the "raw modal data" (quadrature component of resonant 
response) is caused by the presence of off-resonant mode contributions 
in the resonant response. Improvements will not be realized whenever 
the poor orthogonality is due to having a poor mass matrix or to an 
inadequate description of the modes. 

6. Predictions of structural system damping based on experimentally 
determined substructure damping matrices have been verified by direct 


measurement. 
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7. Further development will be required to enable the separation of ' 

closely-spaced modes from the total response at resonance. This is \ 

a problem even when the Imaginary parts of the modes are of no j 

concern; but it is more of a problem when the complex damped modes j 

f 

are being sought because the requirements on frequency separation 
are more stringent. 

8. It may be necessary in applying this method to define modal frequencies 
independently of resonant excitation frequencies because of the 
difficulty in exciting a structure precisely (five or six-digit 
accuracy) at a modal frequency. In this case, additional information 
will be required to help identify the modal frequencies. 

9. Phase errors introduced by the transducer and analog tape systems 
have each been held to + 0.1 8 . If phase errors in the Co/Quad 
data reduction are held to a comparable level (and it appears that 
they are), the total measurement error in the analog system should 
be on the order of + 0.2°. This degree of accuracy should permit 
the measurement of nonproportional damping in systems with damping 
levels as small as 1% to within approximately 30% accuracy. The use 
of digital data acquisition systems should further reduce measurement 
errors in the latter two stages significantly. The basic limitation 
would then lie with the transducers and conditioning equipment. 
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1A. Recommendations for Future Work 

While the basic methodology for damping synthesis by the matrix 
method has been developed and demonstrated, there are certain limitations 
to its current applicability, as already discussed. The following recommenda- 
tions are made for further refinement: 

1. Investigate methods for separating closely-spaced modes. 

2. Investigate methods for more accurately identifying modal frequencies. 

3. Implement the use of digital data acquisition and reduction. 

A. Design and test a physical model without external damping, where 
the substructures' resonant frequencies are more closely spaced. 
Compare predicted and measured values of system damping. 

5. Develop a systematic procedure for estimating the uncertainty 
in predicted values of system damping. This should be based 
on an appropriate statistical characterization of all input 
data. 
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APPENDIX 

CONVERGENCE OF THE MODAL SEPARATION PROCEDURE 


The identification of damped modes depends on being able to implement 
an iterative procedure for removing the off-resonant mode reponse from the total 
response of a structure at resonance. In order for the process to converge, cer- 
tain conditions should be satisfied. While an exhaustive investigation of this 
matter has not been attempted, some relationships have been derived which may be 
interpreted as criteria for convergence. They are not claimed to constitute suf- 
ficient conditions or even necessary conditions in a strict sense. They are 
plausible, however, and seem to be consistent with experience. 

In deriving tnese relationships, it will again be convenient (and 
will avoid confusion) to consider the response vector x to represent flexible-body 
motion only. Thus, when discussing the total response of a free-free structure 
excited by forces which do not add vectorially to zero, it will be assumed that 
rigid-body response has been removed. 

In the present study, vibration tests were conducted on "free-free" 
structures which were externally damped. Strictly speaking, they are no longer 
free-frue in this case. Coupling betweer the rigid-body modes and flexible-body 
modes is introduced by the dampers. Thus the matrix c*^®, which appears in 
Equation (19b) of Section 6, has nonzero elements. This need not complicate the 
study of convergence. Although forces arising from this term may be included in 
the generalized forces which excite flexible modes, it is simpler, and does not 
sacrifice generality to ignore them altogether. Besides, this type of damping is 
artificial and will presumably not be encountered in most practical applications. 

Perturbation Analysis 

It is assumed that the n equations of motion for a real structure may 
be written in the form [Equation (31) of Section 7] 

px(t) + px(t) + <x(t) * f(t). (Al) 

Consistent with the approach taken in [6], the n second-order differential equa- 
tions of (Al) may be put in first-order form leading to a 2n elgenproblem. The 
solution of this elgenproblem will result in a set of complex eigenvalues X and 
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eigenvectors , where X is considered to be a complex diagonal matrix, and a 
complex full matrix. The elements of X are of the form 


' + l "i 

where is the modal decay rate for mode j, and is the damped natural fre- 
quency of that mode. The jth complex eigenvector may be expressed in the form 

♦ . - + i6 *T (A3) 

j R j R j l i 


where $ is defined to be the classical undamped mode. The equations of motion 
R j 

may be transformed according to 


♦r 11 


<A4> 


where 


a is a real matrix whose columns are d> 
R 


Then, (Al) becomes 


Iq + Cq + a) Q q 


♦r £ <‘> 


(A5) 


where it is assumed that <f> R has been normalized so as to result in <p R y I 

and dJ k d> ■ Here 1 denotes an identity matrix while u? denotes a real diag 

R t R o o 

onal matrix whose elements correspond to the undamped frequencies squared of the 

T 

various modes. The matrix £ * P a full modal damping matrix, i.e., non- 

R R 

diagonal in general. It will be assumed that f(t) is of the form 


f(t) - P v g,(t) 
j 3 


(A6) 


where P*. is a vector characterizing the spacial distribution of the forcing func 
J 

tion f(t), used to excite the jth mode, and g(t) is the corresponding scalar 
function of time. Then, (A5) may be written 


Iq + £q + o> o q 


P n MO 
"j i 


(A7) 


where P 


- p . 

’j R x j 



■**w\*:- 
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It will be convenient to represent £ as the ..urn of two matrices 

C - £ + £ (A8) 

a n 

where is a diagonal matrix containing the diagonal elements of £ , and is 

a matrix with a null diagonal containing all the clf-diagonal elements of 
Thus, may be considered to be the diagonal damping matrix of the system. Equa- 
tion (A7) then becomes 

Iq + (£ d + 5 n ) q + <*> 2 q - gj(t) (A9) 


The Laplace transformation of (A9) leads to 


[(ifi) 2 I + (ifi)(£ d + £ n ) + u> 2 ] H q ^ (in) - 


(A10) 


where H is the frequency response vector in the q coordinate system. Further 

q 

defining a complex Impedance matrix Z(ift) by 


Z(ifi) - Z.(ifi) + Z (in), 
a n 


(Alla) 


z d (m) - [(u> 2 - n 2 ) + (in) £ d ] , 


(Allb) 


z (in) - (in) £ , 

n n 


one may express (A10) in the form 


[z.(in) + z Un)J H (in) - p„ 
d " ’l "j 


(A12) 


Since Z. (m) is a diagonal matrix, a scaling transformation 
d 

q - Z~ l/2 (in)Y 

may be made leading from (A12) to 


[I + Z ] H - P , 

n Yj 


(M3) 


where Z - Z~ 1/2 Z Z~ 1/2 , H - zj' 2 H and P - Z~ lf2 P . Then 
n d n d Y J d q j Y J d q j 


H - [I + Z 1 P 
Y 4 « Y< 


(A14) 
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Examination of Z n reveals that lta diagonal elements will be zero, and that its 
off-diagonal elements will be much less than unity in magnitude whenever 
Q*pj ■ u 0 for any frequency u>oj » provided that adequate separation exlst9 among 
the varlou^ reBonant frequencies Of the structure. For these discrete frequen- 
cies, (A14) becomes \ 


H (i P.) - [I + Z (i P.)] P (i P.) 
Yj J i 1 Yj J 

The matrix [I + Z (i P )] * has the series representation 

n J 


(A15) \ 


[I + Z 


n («»J ]" 1 - I-z n (i<V 

n J k-o n J 


(A16) 


provided that the eigenvalues cf Z n (iP^) are less than unity in magnitude [10]. 
In practice, this condition may be established on the basis ' f Gershgorin's disk 
theorem [11] which states that all the eigenvalues of the complex matri.. G lie 


in at least one of the disks of radius r 


j " >: 

J Ul 


G.. centered at G... It is 

1 Jk 1 JJ 


clear that on the basis of (A16), (A15) may be expressed in the manner 


H y (ifij) - [I - Z n (ifJ j )] P y (i^) 


(A17) 


to first-order approximation. Consistent with this formulation, it is usoful to 

also expand H (ip. ) and P (in.) in power series about some nominal values 
Yj J Yj J 

H (in.) and P (10), retaining only the first two terms in each case. Equa- 

Y °J Y °J 

rion (A17) then becomes 


[h (m.) + 6H (in.)] - I - z (in ,)| f 

[\ J ’1 J j [ ” J j[ 


P (ip.) +6 P 
Y Oj J Y J 


(A18) 



.it**- 
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Treating H , 1, and P as zeroth-order terms and 6H , Z , and <5p as 
Y °J r tj y } n Tj 

first-order tert»s, one may equate terms of the same order resulting In 


a (la.) - p v (m.) 

°J °1 


6H (in ) - 6P (ia ) -z (in ) p (in ) 

Y j 3 y ) 3 3 y o j 3 


Transformation from the y back to the q coordinate system, with the notational 
convention Z(i£L) = Zy gives 


H - Z" 1 P 
q d, Q 

°3 j \ 


(A19a) 


6H * Z" 1 6P - Z" 1 Z Z" 1 P 

q j d j q i d j d J q c 


(Al9b) 


In order to define P a , it is required that some initial approxima- 

°3 

tlon of the undamped modes be available. The quadrature acceleration response 
provides this information. The normalized acceleration response corresponding 
to the total response of the structure whe^ excited by a sinusoidal force at fre- 
quency is therefore considered to be the sum of two components 


where 


Hj. * ~ j^CO(xj) + i QUAD(xj)J 


g(t) - g Q sin 8jt 

x j •* x Q sin (fijt + e^) 


CO(xj) * x Q cos 0^ 


OUAD(x, ) ■ -Q, x sin 9. 
J j o j 


Initially, then, ^ may be assigned the values 

L - QUADCH” )/|QUAD(H; ) T y QUAD(H” )J 1/2 

Rj Xj Xj Xj 


(A20) 


«wK4iTk«nnp 



1 




t ! 

•» *w. ”»m r ;...w, ^ v .+ .- w -. - , ■ -v-a 




■ II ." "" - 
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Then, P a is defined to be 

q °j 

P - 2 P 
q 0j R X j 

Since Z, (in ) has been defined in (Allb), H (ifl.) is also now defined. 

d J °j J 

The vector P has bean used in (A7) to represent the total general- 

q j 

ized force vector in the q coordinate system. With P considered to be known, 

T ^ 

P ■ <$> P is determined by <J> D . Then 

q .i R x j R 


S %- % - */ ^ 

The interpretation of 6P as a first-order term in (A18) is valid as long as 


’j 




I ,-1/2 ,T 
^ * 


« 1 


(A21) 


j 


An alternative form for Equation (35) of Section 7 applied to the 
first iteration of this procedure, is 


* R - *(H D) 
R q o 


-1 


(A22) 


where <p * H x D and D is a complex diagonal matrix having elements 


D.. - (e T . H e.)" 1 

jj j q Q j 


(A23) 


The vector e, corresponds to the jth column of the identity matrix. The matrix D 

J * 

acts as a normalizing matrix such that wh^n $ represents response in only the 

A A J 

resonant mode, then - R(<p . ) * $ to first-o. der approximation. From (A22), 

R j J R j 


R(<t> - *) * RUtD" 1 H" 1 -I)]. 

q o 


(A24) 


\ 


I 




- 4 
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o 


Since the diagonal elements of the matrix H D are, by definition, unity, one 

q o 

can represent H D as 
q o 


H D - I + A 


(A25) 


where A is a matrix whose diagonal elements are zero. If the eigenvalues of A 
are less than unity in magnitude, a series representation may be used for 
(H_ D) _1 as discussed earlier. Furthermore, with adequate frequency separation 
of the modes, the elements of A will be small in magnitude compared to unity. 

In this case, the ipproximation 


(H D)” 1 = I - A 
q o 


(A26) 


is valid. Substitution of this result into (A24) leads to 


R($ - *) - Rt«Kl - A - I)] = R(-$A) 


(A27) 


Use of (A22), (A26), and (A27) in (A21) gives 


| 2 - 1/2 R(-»6) T P x 


X 


1 


l in RU(I-A)] 1 P | 

j j 


< < 1 


(A28) 


It is clear that (A28) will be true if the elements of A are in fact small compared 
to unity. 

The element A. . of A is 
kj 


■ \ <v> v 


(a;9) 



1 


1 


1 


<m 4 
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From (23), it is found that the jth element of D Is 


\ v l ■ <•; \ )_1 

i 


2 d l <“i> ♦r P x j >' 1 


A T 

V V 


<% } - n )> + “ 0j > 


When £Jj * u o , then 


1 < 2 0 “o 

1 


"T 

♦r P x 

j j 


The vector (H D)e. * I).. H 1 b found to be 
q c J JJ q 0j 


D.. H 

a <», 


d jj z d 1<iQ J » *i \ 


Because (iQ^ ) is a diagonal matrix 




Dr, H 


1 < 2 ^> 

J 


£ p 

j x j 


for ft. * u> . Then 

3 °J 


ua c ) -j 

w 

O, 


*T 

V»« 

( \ - + 1(2t k \ “ 0j >’ 


(“<) •* )]V-, 


25 J± ) JJ P x 
V k \ / 3 3 



» 




' ?5+Yc»V * -»? t 




--J. 
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It Is clear from this result that to keep the magnitude of A^j small compared to 
unity, it is required that 




172 


< < 1 


(A32a) 


It is meaningful to denote the quantity 




(A3 2b) 


where Q could be interpreted as a normalized dynamic amplification 

facto f->r the kth mode befng excited at frequency w . Thus, (A32a) may be writ 

°j 

ten in the form 




< < 1 


(A32c) 


where P 
excites the 


represents the kth element of the generalized force vector which 
structure at frequency w 0 * The requirement will be satisfied if 




i 


1 
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there Is adequate frequency separation. Otherwise, it will be necessary to shape 
the force distribution P so as to achieve 

X j 


♦: p 


A x i 

— 1 < < i 


♦ » p 


(A33) 


^ *j 


The justification for ignoring in the first iteration depends on 
being able to write 


U + z n (ity] : (i - z n (i^)] 


(A34) 


It is therefore required that each element e^ Z R e^ of Z^ satisfy the relation 


e, Z e. < < 1 

1 k n j 1 


(A35) 


From the definition of Z which follows Equation (A13) , 

tl 

z - z " 1/2 z z " 1/2 

n d n d 


2 n e k 4 ^ .J ^ 


(A36) 


It is recalled that 


T 7~ l/2 e 
e j Z d ®j 


Km 2 - a 2 ) + i2c, « ol “ 1/2 

°j j j 


and that 


T 

ei Z e, 
k n j 


(A3 7) 


where g is the kjth element of the modal damping matrix g for j^k. The largest 
kj 





> * 
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value of e, Z e. is realized whenever ft equals either u or u Without loss 
k n j ° k °j 

of generality, it nay be assumed that ft ■ to < to . Then 

°J °k 


T -1/2 ^ _ ,,, r 2 .-1/2 _ f . 0r 2 .-1 ... r 2.1/2 

e, Z A e j ■ (12 c. to ) - (12c . w ) (i2c to ) 

j j j j «j i i »j 


Furthermore, 


i « *L (i2c, ) 1 

°j °kj J °J 


C, /2C u 
kj j °j 


(A38) 


{ ! 


Substitution of (A37) and (A38) into (A36) for ft ■ to leads to 

o. 


t : 

e, z e J 
k n j 


i2 ^ “o 


1 


2 2 

(to - to ) + 12?, to to 

°k °j k °J °k 


1/2 


kj 


2 C. u 


3 °4 


1(2^) 


2 

o. 


i - 




u> 


1 1 2 ^k - 1 

1 k (i> 

°k 


1 - 


2 

CO 

2 

10 

°k 


1/2 


(’•- !)’ 


n, 




2t. » 

J 


? 'f 


to 


2C 


i t ; 

[ e, Z e. 

1 k n j 1 


j 2 

to 

O, 


1 - 


iV2 


2.2 / \2 

<) *(’“\) 


>72/ 




2 b “o. 


(A39) 
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0 


The use of Equation (28a) from Section 6, and the use of (A32b) in (A39) suggest 
that (A35) may be written 


•ji-iV 172 


iy. 

■jj 


< < i 


(A40) 


With regard to frequency separation, Equation (A40) is seen to be more restrictive 
than (A32c) because is, for the most part, expected to be less than unity. 

A tentative criterion for frequency separation can easily be derived 
by writing (A40) in the form 

1/2 

<Q kj ) - B < < 1 

where B - B|?../5. . | and B - |e. 2 e . | . Squaring both sides of the equation and 
J j k j K n j 

defining to /to - 8 > 1, gives 
°k °j 

4 ?^ 2 2 2 2 

(6 Z -D^ + 4 C 


B 4 

B 4 - 2(1 - 4c *> g 2 + 




For small c^. 


B 2 2 


1 + ^-( l -T*) 




If B is interpreted as the largest value permissible for any |e fc e |, all k,j, 

then for to /to > 1, it would be required that 
°k °j 

n 


^ f¥W) 


: all j ,k 


(A41) 
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For example, not knowing ahead of time, one might tentatively assume 

that £ 1* Then It would be required that 




*cr 




For c ■ .01, B ■ .2, the minimum ratio u /w permitted would b?. 

°k °J 


0) 


— - \[l + ■=# ' 1.225 

w W -04 

°j * 


However, the value of B * .2 was selected somewhat arbitrarily, and will, in fact, 
depend upon both the size and distribution of the scaled impedance matrix Z. The 
larger the dimension of Z (i.e., the larger the number of modes considered), the 
lower one would expect the upper bound B to be. Another way to look at (A41) is 

rp *■ 

to express I ej; Z e. | as a function of (?. 4 /5..), 6 and ? . 

k n j k j j j j 


|e* Z e 
k n j 



(A42) 


Thus, the largest absolute value of any element of the scaled impedance is se?n 
to be 


Z e 1 


fisi 


n j 'max 

; 

hi 

\ (B 2 - 1) 

m 


max 


\e! Z e | 

' k n j 'max 


must be in order to ensure convergence 


The question of how small 

has not been investigated. This condition appears to be fairly restrictive for 

many conceivable applications. However, it could turn out that modes which are 

closely spaced in the frequency spectrum have very weak coupling in the damping 

matrix. If two modes are orthogonal mainly because they represent local motion 

in different parts of the same structure, one would expect corresponding values 

of to be small. Another example is given by a beam which has weak coupling 
kj 

between in-plane and out-of-plane vibration, but nearly the same frequency in 
each case. The isa^lications nere suggest a need for further investigation. 



, 1 . J . : 
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So far in this treatment, it has been assumed that the values of £ 
are deterministic. In practice, this will not be the case since there will always 
be some element of uncertainty. A statistical treatment of the problem has not 
yet been attempted. One would expect that the implication of randomness in £ 
might have an impact on the problem of convergence. If so, the requirements 
involving exciter placement, frequency separation, amount of damping, and measure- 
ment accuracy will be interrelated. A systematic accounting of randomness may 
be of value at some future time, after other refinements have been made, and more 
experience with the method is acquired. 



